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Abstract 


The  data  base  previously  developed  for  multicomponent 
SIALON  Ceramic  phase  diagrams  has  been  expanded  to  cover  Ce90 
additions.  The  Ce203-Si02  and  Ce^-Si^  systems  have  been  3 
computed.  Isothermal  sections  in  the  MgO-Si^N^-SiO  , 

Y2°3“Sl02~Si3N4  and  Ce2°3~Si02"Si3N4  systems  near  2000°K  have 
been  computed  and  compared  with  limited  experimental  data. 

The  trajectory  of  ordering  temperatures  for  A2/B2  and  B2/D0 
reactions  has  been  computed  along  the  Fe3Si-Fe3Al  composition 
path  in  the  BCC  of  the  Fe-Al-Si  system  and  compared  with 
experiment.  The  two  phase  (fee  +  bcc)  fields  for  disordered 
phases  in  the  iron-aluminum-nickel,  iron-aluminum-manganese, 
and  the  iron-nickel-manganese  system  between  900  and  1200°C. 
Construction  of  a  data  base  covering  fluoride  systems  consist¬ 
ing  of  systems  containing  ZrF^  which  are  employed  to  synthesize 
fluoride  glasses  has  been  initiated. 
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I.  PROGRESS  DURING  THE  CURRENT  YEAR 

1.  Calculation  of  Ceramic  Materials  Systems  (Work  statement 
items  g,h,i,j,k)  A  data  base  is  being  developed  for  calculation 
of  quasi-binary  and  quasi-ternary  phase  diagrams  of  ceramic 
systems  (1-3; .  Previous  segments  of  this  base  cover  combina¬ 
tions  of  Cr203,  KgO ,  Al2°3'  Fe2°3'  Fe3°4'  "Fe0" »  Si02,  CaO, 

Si^N^  and  AIN.  Lattice  Stability,  Solution  and  Compound  Phase 
Parameters  were  derived  covering  the  liquid,  spinel,  corundum, 
periclase,  crystobalite ,  tridymite,  quartz,  hexagonal  and  beta 
prime  phases  which  appear  in  the  binary  systems  composed  of 
pairs  of  these  compounds.  Compound  phases  formed  from  specific 
binary  combinations  of  these  compounds  (i.e.  MgO*Cr2C>3)  were 
also  characterized.  This  description  is  based  on  observed 
thermochemistry  and  phase  diagrams  for  the  binary  systems  of 
interest.  Selected  ternary  systems  have  been  computed  based 

on  the  foregoing  data  base  for  comparison  with  experimental 
sections  in  order  to  illustrate  the  usefulness  of  the  data  base. 

The  present  paper  extends  the  data  base  to  cover  BeC,  Y2C>3  and 
Ce203  additions.  Moreover,  ternary  sections  in  the  SiC^-MgO-Si^N^ , 
Si02-Y203~Si3N4  and  SiC^-Ce^-^-Si-jN^  were  calculated  between 
1900K  and  2100K  for  comparison  with  experiment.  Pages  5-26 
comprise  a  paper  published  in  CALPHAD  5^  ( 1 9 8 1  >  185  covering 
this  work. 

2 .  Calculation  of  1II-V  and  II-VI  Electro-Optical  Materials 
Systems  (Work  statement  item  e)  A  data  base  covering  the  binary 
systems  composed  of  Aluminum,  Gallium  Indium,  Phosphorus,  Arsenic 
and  Antimony  has  been  constructed  by  analyzing  the  fifteen 
combinations  of  these  elements  in  terms  of  lattice  stability, 
solution  phase  and  compound  parameters.  Partial  isothermal  sec¬ 
tions  in  the  P-In-As,  As-In-Sb,  P-Ga-As,  Ga-Sb-In  and  Al-Sb-Ga 
systems  were  then  calculated  using  the  foregoing  data  base  * 
comparison  with  experimental  isothermal  sections  and  quasi-binary 
III-V  phase  diagrams.  It  was  found  that  ternary  liquid  and 
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III-V  compound  interaction  parameters  were  required  to  attain 
'  good  agreement  in  some  cases.  Similar  calculations  were  per¬ 

formed  for  the  Te-Cd,  Hg-Cd  and  Te-Hg  binary  systems  and  the 
Cd-Te-Hg  ternary  systems  at  pressures  up  to  74  atmospheres. 
Comparison  of  the  calculated  results  with  experimental  data 
1  on  tie-line  compositions  between  Cd-Te-Hg  liquid  and  quasi¬ 

binary  CdTe-HgTe  alloys  is  important  in  the  liquid  phase  epi¬ 
taxial  growth  of  controlled  band  gap  electro-optical  materials. 
Pages  27-57  comprise  a  paper  published  in  CALPHAD  5  (1981) 

*  217  Covering  this  work. 

3 .  Calculation  of  Ordering  Temperatures  for  the  Fe ^Al-Fe ^Si 
Quasi-Binary  Section  (Work  statement  item  o)  Inden's  model 
has  been  employed  to  compute  the  trajectory  of  ordering 
temperatures  along  the  quasi-binary  join  Fe^Al-Fe^Si  for  A2/B2 
and  B2/D0^  reactions  in  the  bcc  phase  of  the  Fe-Al-Si  system. 

The  results  are  compared  with  experimental  findings  on  page 
67,  while  pages  58-71  detail  the  method  of  computation. 

i 

4 .  Calculation  of  the  Two  Phase  (FCC+BCC)  Fields  in  Disordered 
Fe-Al-Ni ,  Fe-Al-Mn  and  Fc-Ni-Mn  Alloys  Between  900°C  and  1200°C 
(Work  statement  item  p)  The  ManLabs  Data  Base  has  been  employed 

>  to  carry  out  these  calculations  which  are  shown  on  pages  73-75 

and  detailed  on  pages  72-83. 

5 .  Development  of  a  Data  Base  for  Calculating  the  Phase 
Diagrams  of  ZrF^  Systems  used  in  the  Synthesis  of  Fluoride 

*  Glasses  (Work  statement  item  1)  Coupled  thermochemical  analyses 
of  the  ZrF4~BaF2,  ZrF^-LaF. ,  ZrF4~PbF2 ,  ZrF^-NaF,  ZrF^-KF, 
ZrF^-RbF  and  ZrF-CsF  binary  system  have  been  carried  out. 

These  results  can  be  employed  to  locate  the  composition  of 

*  maximum  liquid  stability  for  synthesis  of  fluoride  glasses. 

H.  PROGRAM  PERSONNEL 

Technical  activities  during  the  past  year  have  been 
}  carried  out  by  L.  Kaufman,  F.  Hayes,  E.P.  Warekois,  J.  Kozaczka, 

J.  Nell,  K.  Taylor,  D.  Birney,  C.  Biddell  and  S.  Drake  of 
ManLabs,  Inc. 
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•m.  TECHNICAL  LECTURES 

Technical  lectures  covering  the  work  of  this  program 
were  presented  at  the  Washington,  D.C.  meeting  of  the  American 
Ceramic  Society,  May  1981,  CALPI1AD  X,  Vienna  Austria  July  1981 
and  the  Conference  on  the  Chemistry  of  Materials  at  High 
Temperatures,  AERE,  Harwell,  England  September  1981.  Dis¬ 
cussions  covering  the  cu-rent  work  were  held  with  personnel 
from  ESM-RADC  Hanscom  Field,  AFOSR ,  WPAFB,  UCLA  Los  Angeles, 
California  and  the  Max  Plank  Institute-Stuttgart  and  Dusseldorf. 
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CALCULATION  0!;  QUASI  BINARY  AND  QUAS I  TERNARY 
OXYNTIRIDF  SYSTEMS- IV* 


Larry  Kaufman,  Fred  Mayes  and  Dunbar  Birnie 
ManLabs,  Inc.,  21  brie  Street 
Cambridge,  Massachusetts  02139  USA 

(This  paper  was  presented  at  CALPMAD  X,  Vienna,  Austria,  July  1981) 

ABSTRACT.  A  data  base  is  being  developed  for  calculation  of  quasi-binary  and  quasi-ternary 
phase  diagrams  of  ceramic  systems  (1-3).  Previous  segments  of  this  base  cover 
combinations  of  Ci^O,,  MgO,  A170,,  Fe.,0^,  l:c,0^»  "PeO",  Si07,  CaO,  Si^N^  and' AIN. 

Lattice  Stability,  Solution  arid  Compound  Phase  Parameters  were  derived  covering 
the  liquid,  spinel , corundum,  periclase,  cr\ stobal i te ,  tridymite,  quartz,  hexagonal 
and  beta  prime  phases  which  appear  in  the  binary  systems  composed  of  pairs  of 
these  compounds.  Compound  phases  formed  from  specific  binary  combinations  of 
these  compounds  (i.e.  MgO,Cr7OJ  were  also  characterized.  This  description  is 
based  on  observed  thermochemistry  and  phase  diagrams  for  the  binary  systems  of 
interest.  Selected  ternary  systems  have  been  computed  based  on  the  foregoing  data 
base  for  comparison  with  experimental  sections  in  order  to  illustrate  the  usefulness 
of  the  data  base.  The  present  paper  extends  the  data  base  to  cover  BeO,  Y70,  and 

Co70,  additions.  Moreover,  ternary  sections  in  the  Si07-Mg0-Si  N  ,  Si07-Y7(Y,-Si.N4 

and  Si07-Ce70^-Si _N j  were  calculated  between  1900K  and  2100K  for  comparison  with 

experiment . 

1 .  J_n Product  ion 

>  'I he  utility  of  computer  bused  methods  for  coupling  phase  diagrams  and  thennochemical 

c(ai-'  for  metallic  systems  has  been  well  documented  in  many  papers  published  in  this  journal. 

considerable  effort  is  being  applied  toward  developing  an  extensive  base  for  metallic 
cisterns.  Recently  similar  efforts  have  begun  in  order  to  provide  a  similar  facility  with 
ceramic  systems.  The  expanded  studies  of  SIALON  composites,  combining  silicon  and  aluminum 
t rides  with  oxides  of  silicon,  aluminum,  magnesium,  beryllium,  cerium  and  yttrium  and 
>{, flier  metals  has  provided  additional  motivation  for  predicting  multicomponent  phase  diagrams 
cjr ceramic  systems . 

2.  Description  of  the  'Ihcnmnhemical  System 
Pmployed  to  Glia ract prize  So 1  tit  ion  and  Compound  Phases 

I  The  method  utilized  fc  describing  solution  and  compound  phases  is  the  same  as  that 

<?/>/!> loyed  earlier  (1-3)  comprising  some  symbolic  usage  which  facilitates  data  handling  as 
indicated  below.  The  free  energy,  (.',  of  a  liquid  (solution)  phase,  L,  in  the  binary  system 
?o.^-BeO  is  given  by  liquation  (1)  where  T  is  in  Kelvins,  R=S.314  J/g.  at°K,  while 

is  the  atomic  fraction  of  BO(i.e.  BO=QBeO)  and(l-x)is  the  fraction  of  YO(i.e.  1  / SY^O-. ) . 

mass  basis  is  thus  one  mole  of  atoms  (i.e.  a  gram-atom) 

This  work  has  been  sponsored  by  the  Air  Force  Office  of  Scientific  Research,  Bolling  AFB, 
Washington,  D.C.  under  Contract  F44620-SO-C-0020 

On  leave  from  Department  of  Metallurgy,  University  of  Manchester-UMIST,  Manchester,  England. 
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(1) 


t.  ^<.tc0Gyq  antf'°GgQ  ai’e  the  free  aiergies  of  one  gram  atom  of  pure  liquid  YOandBO.  Table  1  defines  ; 

-//)e  lattice  stabilities  of  the  liquid  and  solid  forms  of  YO,  BO  and  the  other  components  of 
Current  interest.  These  data  derived  from  earlier  studies  and  compilations  of  thermochemical  j 
^Wd  phase_ diagram  data,  (1-8),  when  combined  with  the  solution  and  compound  phase  parameters 
vS  town  in  Tables  2  and  3  permit  calculation  of  the  binary  phase  diagrams  shown  in  Figures  1-12. 
)4hc  solution  parameters  LYOBO  and  l.BOYO,  which  describe  the  liquid  Y.O.-BeO  solution  are 
listed  in  Table  2  (LYOBO=LBOYO=-2092  Joules/g.at) .  Similar  parameters  ror  the  liquid  and 
'S’Olid  phases  are  listed  in  Table  2.  The  free  energy  of  the  solid  phases  are  described  in  a 
/4/<?nner  similar  to  equation  (1).  Thus,  for  example,  the  free  energy  of  the  body  centered 
Cubic  (MUjO^JY  phase  in  the  YO-BO  system  is  given  by 

>  <-=(l-x)°GY0+x<bp0+RT(x£nx+(l-x)£n(l-x))+x(l-x)  [YYOBO(l-x)+xYBOYO]  J/g.at.  (2) 

•/•/? e  solution  parameters  for  the  Y  phase  (YYOBO=YBOYO=-2G92  J/g.at)  are  listed  in  Table  2.  The  free 
energy  difference  between  the  I,  and  Y  forms  of  YO  and  BO,  i.e.  YOYOLY=°G|;q-0G^q  and 

4-°GYb0=B0B0LY  are  listed  in  Table  1.  The  free  energy  of  a  compound  phase  such  as 

*  =(1/7)  (SiOo*2Be0)=SO0  571  *s  tlol‘"inCc*  on  Basis  ot"  iaBle  3  in  terms  of  the  base 

Masc ,  the  compound  parameter,  C,  and  the  stoichiometry  by  equation  (3) 

&=0.429Ggo+0.S71GgO+(0.429) (0.571) [0.429LSOBO+0.571LBOSO  -C]  J/g.at.  (3) 

ife-forence  to  Table  2  shows  that  LSOBO=bBOSO=41840  J/g.at  while  Table  3  shows  that  0=62132-7. 406T. 

y/jus 


^»0.429GgO*0.571Gjo-4971*1.813T  J/g.at. 


(4) 


At  29SK  reference  to  Figure  1  shows  that  the  stable  form  of  quartz  is  the  beta 


iform  (designated  R) .  Table  1  shows  that 
50'GS0=S°S0,‘R'S0S0LN=5042‘1  ’ 0041  J/g.  at , 


(S) 


substitution  of  equation  (5)  into  equation  (4)  yields  an  expression  for  the  free  energy  of 
formation  of  (1/7) (SiO  • BoO)  at  298K  from  the  stable  forms  of  it  s  compound  components  as, 

S,.l298Kj=GG-0.429Gg0-0.571G|Jn=-2808+1.382T  J/g.at. 


BO 


(6) 


Reference  to  Table  4  shows  that  this  result  is  in  keeping  with  the  experimental  (8) 
■/fllermochcmical  data  on  the  free  energy  of  formation  of  2BeO*SiO,,.  Table  4  also  displays  the 

Iculated  free  energy  of  formation  for  the  remaining  compounds  listed  in  Table  3  and  shown  in 
"§gures  1-12.  It  should  be  noted  that  this  compound  parameter  for  SNn  „SO  currently  shown 

AS  0=62760-1 .2S5T  J/g.at.  has  been  revised  from  the  previously  stated  value  C=1 15060-25. 10T 
/g.at.(3)fhe  current  value  is  based  on  the  assessment  of  Doe^ner  et  al  (5)  and  leads  to  slight 
evisions  in  the  previously  calculated  SN-SO  and  SO-AO-SN  phase  diagrams  presented  earlier. (3). 
revised  versions,  which  differ  slightly  from  the  earlier  results  are  shown  in  Figures  8 
13. 

3.  Calculation  of  Quasi-Ternary  Phase  Diagrams 

The  free  energy  of  ternary  solution  phases  was  synthesized  from  the  binary  solution 
phases  on  the  basis  of  Kohler's  equation  as  in  the  previous  papers  (1-3).  On  this  basis  the 
kfroe  energy  of  the  liquid  phase  in  the  SO-MO-SN  system  where  x  is  the  atom  fraction  MO,  y  is 
fhe  atom  fraction  SN  and  1-x-y  is  the  atom  fraction  SO  is  given  by 
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TAB  LI-  1 


SUMMARY  OF  L AIT  ICE  STABILITY  PARAMETERS 
(All  units  in  Joules  per  gram  atom  (mole  of  atoms),  T  in  Kelvins) 


AO 

ft 

BO 
P  = 

R  a 

I 


=  (1/5)A1,03,  SO  =  (!/3)Si02,  MO  =  (l/2)MgO,  AN  =  (1/2)A1N,  SN  =  (l/7)SijN4 
•  (l/2)BeO,  YO  =  (1/5)Y203,  CE  »  (l/S)Ce203 

Periclase,  C  =  Corundum,  X  =  Crystobalite,  T  =  Tridymite,  H  =  a  quartz 

C  quartz,  B  =  beta  SijNj,  N  =  hexagonal  AIN  and  BeO,  Y  =  body  centered  cubic  (Mn203)Y203 
and  Ce203  structure 

BOBOLN*  =  (1/2) BeO( Liquid) - (1/2) BeO(hcxagonal) 

YOYOLY  =  ( 1/5) Y-0_( Liquid) - (1/5) Y^O,  (body  centered  cubic) 

*  w  C  J 


ft 


YOYOLY 

3 

22694 

- 

8.368T 

C  EC  ELY 

= 

20334 

- 

8.368T 

Y0Y0LN 

S 

0 

- 

14.142T 

CECELN 

s 

0 

- 

14.142T 

YOYOLX 

= 

0 

- 

2.092T 

CECF.1.X 

= 

0 

- 

2.092T 

Y0Y0LB 

0 

- 

12.51 0T 

CECELB 

- 

0 

- 

12.51 0T 

YOYOLC 

= 

0 

- 

10.209T 

CECELC 

= 

0 

- 

10.209T 

BOBOLN 

3 

40376 

- 

14.142T 

A0A0LY 

0 

* 

8.368T 

BOBOLX 

ts 

0 

- 

1.674T 

ANANLY 

3 

0 

- 

8.368T 

BOBOLC 

52 

0 

- 

10.209T 

SNSNLB 

= 

33949 

- 

12.510T 

BOBOLB 

3 

0 

- 

12.510T 

SNSNLY 

= 

0 

- 

8.368T 

BOBOLY 

3 

0 

- 

8.368T 

S0S0LY 

= 

0 

- 

2.092T 

SOSOLN 

c 

0 

- 

2.092T 

SOSOLX 

- 

3347 

- 

1.674T 

SOSOLR 

50-12 

- 

3.096T 

*  These  differences  specify  the  free  energy  of  one  phase  (i.e.  liquid)  minus  the  free  energy 
of  the  second  phase  (i.e.  hexagonal)  for  a  given  compound. 


.L  ,,  „  «pL  ...0/.L 


G  =(l-x-y)  Gg0*x  GM0+y°GSN+RT((l-x-y)£n(l-x-y)*x€nx*yfny] 

>(l-x-y)x(l-y)'l[(l-x-y)LS0M0+xLM0S0]*xy(x+y)‘1[xLM0SN+yLSNM0] 

ft  +(l-x-y)y(l-x)“1[(l-x-y) LSOSN+yLSNSO]  J/g.at.  (7) 

The  solution  parameters  required  to  specify  G1"  were  defined  previously  (2,3)  as 
LSOSN=LSNSO=LMOSN=LSNMO= 29288  J/g.at.  and  LS0M0=lM0S0=- 106274+42. 01T  for  0 . 4< x^ 1 . 0  and 

LS0M0=-24267+24.27T;  LMOSO=-229283+68.62T  for  O.O^x^O.IO.  These  parameters  permit  explicit 

^definition  of  GL  in  composition  ranges  where  (x/l-x-y)  is  more  than  0.667  and  ranges  where 
(x/l-x-y)  is  less  than  0.667.  In  the  latter  range  the  miscibility  gap  present  in  the  SO-MO 
binary  system  propagates  into  the  ternary  S0-M0-SN  as  shown  in  Figure  14a.  Here  the  ternary 
miscibility  gap  shows  the  tie  lines  connecting  the  coexisting  liquid  compositions  L^  and  L2» 

The  remaining  solution  phases  in  this  system,  i.e.  P,  X,  B,  etc.  are  defined  in  similar 
g fashion  to  Equation  7.  The  resulting  equilibria  between  the  L  and  P  solutions  are  depicted 
by  the  tie  lines  traversing  the  two  phase  L+P  field  shown  in  Figure  14(b).  The  compound 
phases  F,R  and  BP  which  appear  in  the  SO-MO,  SO-SN  and  MO-SN  binary  systems  are  defined 
along  the  lines  previously  established  (1-3)  and  are  specified  in  Table  3.  Thus  for  example 
the  free  energy  of  the  SNn  _S0n  ,  (R)  phase  in  the  SO-MO-SN  system  is  defined  by  Equation  8 
as  u,/  ' 


TABLE  2 


SUMMARY  OF  SOLUTION  PHASE  PARAMETERS 
(All  units  in  Joules  per  gram  atom  (mole  of  atoms) ,  T  in  Kelvins) 

L  =  Liquid,  N  =  hexagonal,  X  =  Crystobalite,  T  =  Tridymite,  Y  =  body  centered  cubic 

B  ■  beta  Si^,  AO  =  (1/S)A1203>  SO  =  (l/3)Si02,  AN  =  (1/2) AIN,  SN  =  (l/7)Si3N4»  BO  =  (l/2)BeO 

YO  =  (1/S)Y203,  CE  =  (l/S)Ce203 


$LSOBO  »  LBOSO  =  41840 
NSOBO  »  NBOSO  =  83680 
XSOBO  *  XBOSO  ■  83680 
TSOBO  »  TBOSO  =  83680 

LYOBO  a  LBOYO  =  -2092 
&YYOBO  »  YBOYO  *  -2092 
NYOBO  »  NBOYO  =  83680 

LAOBO  a  LBOAO  =  8368 
CAOBO  «  CBOAO  a  83680 
NAOBO  »  NBOAO  a  83680 

41- 

LBOAN  »  LANBO  a  20920 
NBOAN  a  NANBO  «  83680 

LCESN  »  LSNCE  a  4180 
BCESN  »  BSNCE  a  62760 
i.YCESN  »  YSNCE  =  62760 


LBOSN  a  LSNBO  a  20920 
BBOSN  =  BSNBO  =  83680 
NBOSN  =  NSNBO  =  83680 


LYOAO  =  LAOYO  a  12552 
YYOAO  =  YAOYO  =  83680 
CYOAO  =  CAOYA  =  83680 

LYOAN  =  LANYO  =  20920 
YYOAN  »  YANYO  a  83680 
NYOAN  =  NANYO  =  83680 

LSNSO  =  LSOSN  =  29288 
XSNSO  a  XSOSN  =  12SS20 
TSNSO  a  TSOSN  =  1255 20 

LYOSN  »  LSNYQ  =  4180 
BYOSN  =  BSNYO  =  83680 
YYOSN  a  YSNYO  =  83680 


0<xso<0.60* 

LYOSO  =  LSOYO  =  -92048+36. 82T 
0.6<xso<1.0 

LYOSO  *  -216,000+66. 94T 
LS(  YO  =  -9414+16. 74T 
YYCSO  =  YSOYO  a  83680 
XYOSO  =  XSOYO  =  83680 
TYOSO  =  TSOYO  =  83680 

0<xSq<0.60 

LCESO  a  I.SOCE  »  -92048  +  36.82T 
0.6<xS(><].0 

LCESO  +  -216,000+66. 94T 
LSOCE  =  -9414+16. 74T 
YCESO  =  YSOCE  *  83680 
XCESO  =  XSOCE  =  83680 
TCESO  =  TSOCE  a  83680 
*x_~  =  atomic  fraction  SO  in  YO-SO 
&  or  CE-SO 


+RT(y£ny+z£nz-(l-xKn(l-x)!i.AGE  J/g.at. 


(8) 


where  xaXgQ,  y=y^jQ  and  zaz^^.  Since  t.«e  compound  phase  runs  from 


( 


SN0.70S00.30  t0  MO0.70SO0.30’  P*^*"* %)/U-xJ*0  and  x=x;+yp=x*=0.3,  and  z=l-x;-y(l+p)s 
0.7-y.  Since 


AGa  =  (0.3) (0. 7) [O.7LSNSO+0.3LSOSN-C]  =  -7029  +  0.263T  J/g.at.  *  (9) 
and  AGg  =  (0.3) (0.7) [0.7LMOSO+0.3LSOMO-C]  =  -22318  +  8.822T  J/g.at. 

(  from  Tables  3  and  5  and  the  previous  values  of  LSNSO,  LSOSN,  LMOSO  and  LSOMO  (2)  given  in  the 
above  text.  The  excess  free  energy  of  mixing,  AG^,  for  the  compound  is  given  (3)  by 

AGg  =  CAByz/(l-x) 

However,  CAB  is  taken  as  zero  in  all  the  cases  treated  here  thus  Equation  8  can  be  written 
rexpiicitly  as  follows: 

GR  =  (°-70-y)OGsN+0.3°Gg0+yeG®0  -  ( l-(y/0. 7)) (7029-0. 263T) 

-(y/0.7)(22318-8.822T)+RT(y£ny+(0.7-y)£n(0.7-y)-0.7tn0.7)  (11) 

^where  Y-Y^  Since  the  lattice  stabilities  of  SN,  SO  and  MO  (i.e.  SNSNLB,  SOSOLB  and  MOMOLB) 
have  already  been  specified  the  R-L  equilibria  (Figure  14c)  can  be  computed.  Figs.  (14a)  through 


TABLE  3 


Compound 


SUMMARY  OF  COMPOUND  PARAMETERS  FOR  BINARY  SYSTEMS 
(All  units  in  Joules  per  gram  atom  (mole  of  atoms)  T  in  Kelvins) 

Name  Stoichiometry  Stability  Base  Compound  Parameter 


(Joules/g.at.) 


(1/7) (Si02-2Be0) 

G 

S00.429B00.571 

stable 

N 

62132-7. 406T 

(1/17) (3A1203-Be0) 

Q 

A00. 882B°0. 118 

stable 

C 

75312 

(1/7) (Al203-Be0) 

I 

A00.714B00.286 

stable 

C 

76567 

(l/il)(Al20  *3Be0) 

J 

AO0.455BO0.545 

stable 

N 

61923+9. 916T 

(1/1S)(2Y203'A1203) 

U 

Y00.667A00. 333 

stable 

Y 

70291 

(i/io)(-y2o3-ai2o3) 

V 

YO0.500AO0.500 

stable 

C 

48534+20. 083T 

(1/40}(3Y2Os’5A12O3) 

WA 

Y00.375A00.625 

stable 

C 

82425 

(1/5) (Si2N20) 

R 

SN0.7SO0.3 

stable 

B 

62760-1. 2S5T 

(1/12) (Y203*Si3N4) 

S 

Y°0.417SN0,583 

stable 

B 

44769+14. 644T 

(1/8)(Y203-Si02) 

M 

YO0.625SO0.375 

stable 

Y 

-15900+29. 288T 

(l/19)(2Y203-3Si02) 

K 

YO„  e,-S0A 

0.S27  0.473 

stable 

Y 

-30334+36. 8 19T 

(l/ll)(Y203-2Si02) 

0 

YO_  .teSOft 

0.455  0.545 

stable 

Y 

-2218+23. 849T 

(l/12)(Ce203-Si3N4) 

S 

CG0.417SN0.S83 

stable 

B 

44769+14. 644T 

(l/8)(Ce203-Si02) 

M 

Ci:0.625S°0.375 

stable 

Y 

418+24. 058T 

(l/62)(7Ce203-9Si02) 

Z 

Cn0.565S00.435 

stable 

Y 

6276+24. 434T 

(l/ll)(Ce203-2Si02) 

0 

CE0.455S00.545 

stable 

Y 

7531+25. 230T 

0.01(14Si3N4*BeO) 

BP 

SN0.98B00.02 

stable 

B 

83680+4. 184T 

(1/150) (21Si3N4*Si02) 

BP 

SN0.98S00.02 

stable 

B 

71128+4. 184T 

0.0040(35Si3N4-Y203) 

BP 

SN0.98YO0.02 

stable 

B 

■  83680+4. 184T 

0.0040(35Si3N4*Ce203) 

BP 

SN0.98CK0.02 

stable 

B 

83680+4. 180T 

gl4(d)  illustrate  the  steps  which  are  taken  in  sythesizing  the  ternary  system  SO-MO-SN, 

(i.e.  l/3Si07-l/2Mg0-l/7Si_N.)  from  the  component  binary  systems.  Initially  each  of  the 
equilibria  between  pairs  or  solution  and  or  compound  phases  are  calculated  individually.  In 
the  present  case  these  pairs  consist  first  of  the  liquid  miscibility  gap  (L.-L2  I4a) 
followed  by  l-Periclase  in  Figure  14b.  Next  the  L-F  and  L-R  equilibria  which  are  shown  in 
Figures  14a  and  14c  are  calculated.  Overlaying  the  L.-L-,  L-F,L-P,and  L-R  equilibria  shows 

id  ^  ^ 

‘‘‘that  there  are  no  interactions  between  any  of  the  phase  pairs  except  in  the  case  of  the 
L-R/L^L^  interaction.  The  latter,  which  's  shown  in  Figure  14c  results  in  a  three  phase  field 

between  R,  L^  and  L^.  The  boundaries  of  this  three  phase  field  are  the  lines  traversing  the 

Lj*L2  miscibility  gap  illustrated  in  Figure  14a  and  the  tie  lines  which  define  the  R-L 

fcequilibria  shown  in  Figure  14c. 
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1/62 { 7Co_0 


0.565S00.435 


L  (liquid) 


TABLE  4 


CALCULATED  FREE 
(All  Units  in  Joules 


ENERGY  OF  FORMATION 
per  gram  atom  (mole 


OF  COMPOUND  PHASES 
of  atoms),  T  in  Kelvins) 


Compound 

i* 

Name 

Stoichiometry 

Eree  Energy  of  Formation,  AGJ298K] 
From  Component  Compounds  1 

(l/7)(Si02*2Be0) 

G 

SO0.429BO0.571 

-2808  +  1. 382T/- 2809+0. 777T(8) 

(1/17)(3A1  0  *BeO) 

Q 

A00. 882B00. 118 

-2203  -  0.464T 

c(l/7)(Al203.Be0) 

I 

A00.714B00.286 

-2381  -  1.125T 

(1/11) (A1203  *  3BeO) 

J 

A00.455B00.S45 

-2524  -  1.049T 

(1/15)(2Y203*A1203) 

U 

Y00.667A00.333 

-4950  -  0.615T 

l(1/10)(Y2O3’A12O3) 

V 

YOO.S0OAOO.50O 

-2351  -  4.100T 

(1/40)(3Y203-5A1203) 

WA 

Y00.375A00.625 

-7866  +  0.690T 

(1/5) (Si2N20) 

R 

SN0.700S00.300 

-5514  -  0.038T/-5781  +  0.40T(5) 

;;U/12)(Y2<ysi3N4) 

S 

Y00.'«17SN0.583 

-404  -  1.833T 

(l/8)(Y203*Si02) 

M 

Y00.625S00.375 

-15958  +  1.389T 

(l/19)(2Y203*3Si02) 

K 

Y0n  ,07SOft  ... 
0.527  0.473 

-13000  -  0.475T 

(1/11) (Y203-2Si02) 

0 

Y0n  .CPS0A 

0.455  0.545 

-19527  -  2.669T 

(1/12) (Co203-SisN4) 

S 

Cr‘0.417SN,0.583 

-1389  -  1.828T 

(l/8)(Cc203’Si02) 

M 

Cri0.625S00.375 

-19787  +  2.615T 

(l/62)(7Ce203-9Si02) 

Z 

CH0.565S00.435 

-21974  +  2.606T 

(l/U)(Ce203*2Si02) 

0 

CEn  .c-SOA  ... 
0.455  0.545 

-21945  +  2.326T 

'  0.01(14Si3N4-BeO) 

BP 

SN'0.980B00.020 

-423  -  0.113T 

(1/150) (21Si3N4-Si02) 

BP 

SN0.980S00.020 

-720  -  0.102T 

0.004(35Si3N4-Y203) 

BP 

SN0.980Y00.020 

-1104 

0.004(3SSi3N'4*Ce203 

BP 

SN0.980CE0.020 

-1158 

The  calculation  of  the  isothermal  section  at  2’00K  can  be  concluded  by  adding  the  L-BP  equili¬ 
bria  as  shown  in  Figure  14d.  In  this  case,  as  in  the  R-L  case,  there  are  interactions  with 
phase  pairs  which  have  been  previously  considered.  First  there  is  an  interaction  between  L-BP 
and  the  liquid  gap,  Lj-L2>  as  illustrated  in  Figure  14d,  then  there  is  a  second  interaction 

between  the  L-BP  and  L-P  pairs.  Comparison  of  Figures  14b  and  14c  with  14d  shows  that  while 
the  L-R/L-BP  and  the  L-BP/L-P  interaction  pairs  form  stable  throe  phase  fields,  the  L-BP/ 
L^-L^  interaction  pair  which  is  shown  in  Figure  14d  is  metastable.  This  can  be  seen  by 

comparing  14c  and  14d  and  noting  that  the  three  phase  BP+L^+L,  field  is  "covered"  by  the 

R+Lj  +  L.,  and  R+L  fields.  Detailed  calculation  of  the  L-BP  equilibria  can  be  performed  by 

defining  the  free  energy  of  the  BP  phase  on  the  bases  of  Tables  1-4. 
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TABLE  5 
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SUMMARY  OF  COUNTERPIIASE  STOICHIOMETRY  AND  PARAMETERS 
EMPLOYED  IN  TERNARY  CALCULATIONS 


System 

Stable  Phase 
(Name) 

Base 

Counterphase 

Base 

Counterphase 

Parameter 

SO-MO-SN 

S00.429M00.57l'F> 

P 

SN0.429M00.571 

P 

0 

S00.300SN0.700^R) 

B 

S°0. 300M00. 700 

B 

0 

SN0.980S00.020^BP^ 

B 

SO0.SO0MO0.500 

B 

0 

SN0.980M°0.020^BP^ 

B 

M00.500S00.500 

B 

0 

SO-YO-SN 

YO0.625S°0.375f'Mj 

Y 

Y°0.625SNQ. 375 

Y 

0 

Y00.527S00.473(K;| 

Y 

YO0.527SN0.473 

Y 

0 

V00.45SS00.545^ 

Y 

YOQ.455SN0.545 

Y 

0 

S00.300SN0.700^R) 

B 

S00.300'1  °0. 700 

B 

0 

SN0.980Sno.020^BP^ 

B 

S00.500'1  °0.500 

B 

0 

SN0.980^°0.020^BP-) 

B 

Y00.500S00.500 

B 

0 

SO-CE-SN 

CK0.625S00.375W 

Y 

Cl:0.625SN0.37S 

Y 

0 

CR0.565S00.435(Z) 

Y 

CU0.565SN0.435 

Y 

0 

CK0.455S°0.545(O) 

Y 

('E0.455SN0.545 

Y 

0 

S00.300SN0.700^R) 

B 

SO0.300C'r:0.700 

B 

0 

SN0.980SO0.020^BP^ 

B 

SO0.500Cl:0.S00 

B 

0 

SN0.980CE0.020(BP) 

B 

U:0.500S00.500 

B 

0 

Thus  the  equilibrium  between  the  liquid  and  BP  phase  can  be  computed  as  in  Figure  14(d). 
Combination  of  the  pairwise  equilibrium  in  Figure  11  yields  the  computed  isothermal  section 
at  2100K  shown  in  Figure  15.  Isothermal  sections  similarly  derived  at  2000K  and  1900  are  sho*>n 
in  Figure  15  along  with  an  observed  section  due  to  Mueller  (9)  at  1500K.  The  latter  is  in 
keeping  with  the  calculations  shown  at  1900K.  (See  "Note  in  Proof") 

Similar  calculations  were  performed  in  the  SO-YO-SN  and  SO-CE-SN  systems.  The 
results,  which  are  complicated  by  the  existence  of  stable  quasi-ternary  compounds  (listed  in 
Table  6)  are  displayed  in  Figures  15-19.  'Die  SO-YO-SN  and  SO-CE-SN  sections  were  computed  in 
a  manner  similar  to  the  SO-MO-SN  case  in  that  the  systems  were  first  computed  as  if  the  ternary 
phases  were  absent.  Subsequently  the  ternary  phases,  designated  as  C,D  and  E  were  inserted. 

,fn  keeping  with  previous  practice  (1-5)  the  free  energy  of  the  ternary  phases  C,  D  and  E  were 
defined  by  choosing  a  base  phase  and  then  defining  a  compound  parameter  in  conformity  with 
experimental  observation.  In  the  present  case  the  Y  structure  was  chosen  as  the  base  phase 
since  the  ternary  phase  composition  occurs  near  the  Y0  and  CE  corners  of  the  SO-YO-SN  and 
SO-CE-YO  systems  as  is  seen  in  Figures  17  and  18. 

The  free  energy  of  the  ternary  compound  D  in  the  SO-YO-SN  system  Dsfl/lSJY^Si^O^N^3 

!  *°0.100YOl).M>7S!f0.223  H  defi"0d  05 


Compound 


TABLE  6 

DESCRIPTION  OF  QUASI -TERNARY  COMPOUNDS  IN  THE 

Si02"Y2°3‘Si3N4  AND  Si02'Ce2°3'Si3N4  SYSTEMS 

.  '  Free  Energy  of 

Stoichiometry  Compound  Melting  Formation  from 

Parameter  (C)  Point  Component  Compounds 

(Joules/g.at.)  (°K)  (AGf[298Kj) 


C=  (l/5)YSi02N  Y 
!  D=(l/15)Y4Si207N2  Y 
E=(l/21)Y5(Si04)5N  Y 
C=(l/5)CoSi02N  Y 
D»(l/15)Ce4Si207N2  Y 


SO0.150YO0.SO0SN0.350 

SO0.100YO0.667SN0.233 

S00.322Y00.59SSN0.083 

SO0.150Cn0.500SN0.3S0 

S00.100CE0.667SN0.233 

S00.322CK0.595SN0.083 


29096+8. 368T  1953 

13004+16. 736T  2110 

-2155+20. 92T  1996 

31505+8. 368T  1914 

14769+16. 736T  2020 


Joules/g.at.. 


-4 393- 2.41 0T 
-678-4. 127T 


■1 1456-0. 883T 


-5422- 2. 410T 
-1259-4. 127T 


E=(l/21)Ce5(Si04)3N  Y  S00.S22CK0.595SN0.083  3975+20. 92T  2109*  -13585-0. 883T 

*SO0.322CEQ.595SN0.083  dcc°mposos  peritcctical lv  into  Z  and  Liquid  above  2100K 

G  =0*10°Rso*0.6679GyO+0.2338GYN+(0.10) (0.667) (0. 767) _1 [0. 1LSOYG+0.667LYOSO-C] 

1  f°-  ILSOSN+O.  233LSNS0-C] +  (0.667)  (0.233)  (0.90) "1  (0.667LYOSN+ 
u.^LbMO-CJ  J/g.at.  (jg) 

.■>ince  C=  13004 +16. 736T  for  this  phase 

G°=0, 10°CY0+0. 667°GyO+0. 233°Gg^-9094-3.062T  J/g.at.  (16) 

fhTfree  °f  f?rmit*oa  °f  D  from  tho  Y  form  of  SO,  YO  and  SN  is  -9094-3. 062T  J/g.at. 

resofctivoWW  of  formatlon  °f  D  from  the  stable  forms  of  SO,  YO  and  SN  (i.e.  R,  Y  and  B 

1j/gjat.StabUity  ValU6S  8iVe"  in  TablC  >•  With 

G^=0. 10°Gg0+0  667°GY+0.233®G^, -678-4. 127T  J/g.at.  (l7) 


AGfl298K]— 678-4. 127T  J/g.at.  (18) 

VnMn  ,  ThP  free  worjy  and  free  enerS>’  °f  formation  of  the  remaining  compounds  listed  in 
wise' I  ,.  ?,,  !  C(7l,ted  »»  »  fashim.  On  this  basis  the  isothermal  sections  and  pair- 

at  2100  w erS  develoPed-  figure  16  shows  the  computed  sections 

G2  °  20°0K  ln  V*i;\pha;ef  P”5ent  in  the  binar>'  systems  and  the  ternary  compound  D 

?900K  Tn  TTeL  ia)  and.16(b)  3nd  Fi*u*»  17(a)  illustrate  pairwise  interactions  at 
j  J  ilnnv  \ -^eSe  ^™1SC  interactions  combine  to  yield  the  calculated  isothermal 
n  Flanrf  ? \WhlCJ  co,\tains  two  small  liquid  regions.  A  computed  section  at  1823K  given 
.®  .  .  c  - 'o''s  ii°w  these  regions  have  grown  smaller  with  decreasing  temperature.  This 

wM^°?nS  f°d  a2rcement  with  the  observed  section  (10)  at  1823  shown  in  Figure  17(d) 
which  illustrates  the  ternary  compounds  and  the  existence  of  a  small  liquid  field  as  well. 

?inn  ?nnn  isotbcrmal  sections  in  SO-CE-SN  is  displayed  in  a  similar  manm'r 

at  2100,  .000  and  1900k.  At  2100K  the  calculated  isothermal  section  given  in  Figure  18 

^  yIhe.PhTLStab1^  V\ the  cds°  binar>'  systems.  At  2000K  the  calculated  section 
-rnai>  (Ce4Si207N:2)  and  b(C’e3(Si04) ,N)  phases.  Pairwise  equilibria  calculated 
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Calculated  and  Observed  (9)  Isotherrea 
(See  "Note  in  Proof”) 


ated  Isothermal 


2000 


se 


(l/3)SiO,=SO 


m  (l/5)Ce_0  =CE 

3*^4  Z  6 

Calculated  Isothermal  Sections  at  2100°  and  2000°K  and  Pairwise  Equilibria  at 


(l/3)Si0-,=S0 


ts&VIsW1 


at  1900K  areshown  in  Figures  18(a),  18(b),  19(a)  and  19(b).  The  latter  combine  to  yield  the 
ralculated  i'sothermal  section  at  1900K  given  in  Figure  19(c).  The  latter  is  in  relatively 
good  agreement  with  the  observed  subsolidus  equilibria  (11)  determined  between  1923  and  2023K. 

The  observed  section  does  not  show  evidence  of  the  S  or  D  compound  phases  which  was  expected 
but  not  observed  (11). 

‘The  SO-MO-SN,  SO-YO-SN  and  SO-CE-SN  sections  calculated  above  serve  to  illustrate 
ftow  the  synthesis  of  multicomponent  phase  diagrams  can  be  performed  in  support  of  experimental 
Studies  of  practical  problems  which  must  often  be  performed  with  great  difficulty  and  some 
measure  of  uncertainty. 
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Note  in  Proof 

"ote  in  Proof 

H.  L.  Lukas  has  pointed  out  that  the  Gibbs  free  energy  change  for  the  reaction  between 
the  following  solid  phases 

Si3N4  +  4Mg0  -►  Mg2Si04  +  2Mg  SiN2 

is  equal  to  3953  -  8.85T  Joules  or  264  -  0.59  T  Joules/  g.  at.  On  this  basis  the  SN  -  MO  binary 
edge  is  slightly  metastable  with  respect  to  the  formation  of  F  (1/7  (Si02  •  2MgO))and  1/5 
(Mg  Si  N2)  above  446  K.  Thus  the  sections  shown  in  Figure  15  represent  metastablo  equilibria. 


*  CALCULATION  OF  TF.RNARY  SYSTEMS  CONTAINING  III-V 

*  AND  1 1 -VI  COMPOUND  PHASES* 

Larry  Kaufman,  Janine  Nell,  Keith  Taylor  and  Fred  Hayes 
Man Labs,  Inc.,  21  Erie  Street 
Cambridge,  Massachusetts  02139,  USA 

(This  paper  was  presented  at  Calphad  X,  Vienna,  Austria  July  1981) 

ABSTRACT.  A  data  base  covering  the  binary  systems  composed  of  Aluminum,  Gallium,  Indium, 
Phosphorus,  Arsenic  and  Antimony  has  been  constructed  by  analyzing  the  fifteen 
combinations  of  these  elements  in  terms  of  lattice  stability,  solution  phase  and 
compound  parameters.  Partial  isothermal  sections  in  the  P-In-As,  As-In-Sb,  P-Ga-As, 
Ga-Sb-In  and  Al-Sb-Ga  systems  were  then  calculated  using  the  foregoing  data  base 
for  comparison  with  experimental  isothermal  sections  and  quasi-binary  III-V  phase 
diagrams.  It  was  found  that  ternary  liquid  and  III-V  compound  interaction  para¬ 
meters  were  required  to  attain  good  agreement  in  some  cases.  Similar  calculations 
were  performed  for  the  Te-Cd,  Hg-Cd  and  Te-Ilg  binary  systems  and  the  Cd-To-llg 
ternary  systems  at  pressures  up  to  74  atmospheres.  Comparison  of  the  calculated 
results  with  experimental  data  on  tie-line  compositions  between  Cd-Te-Hg  liquid 
and  quasi-binary  CdTc-HgTc  alloys  is  important  in  the  liquid  phase  epitaxial  growth 
of  controlled  band  gap  olcctro-optical  materials. 

1 .  Introduction 

Alloy  semiconductors  arc  important  for  a  wide  range  of  electro-optical  applications 
because  their  properties  can  be  tailored  by  controlling  the  composition  of  the  solid,  which  is 
grown  from  the  liquid  or  deposited  from  the  vapor.  Although  crystal  growers  have  become  very 
adept  at  controlling  the  temperature  (and  pressure)  of  the  parent  phase  in  order  to  obtain 
the  desired  characteristics  of  the  crystalline  alloy  semiconductors,  phase  diagram  data  is 
indispensable  in  order  to  deal  with  the  increasingly  complex  systems.  Stringfellow  and 
co-workers  (1-6)  have  pioneered  in  developing  simple,  useful  models  for  calculating  quasi¬ 
binary  systems  based  on  III-V  compounds,  Stringfel low' s  pioneering  work  has  been  extended 
and  expanded  by  Ansara  and  co-workers  (7),  Osamura  and  Murakami  and  co-workers  ( 8 - 1 0 )y  and 
Brebrick  (11)  who  have  used  more  extensive  experimental  thermochemical  data  and  less  restrict¬ 
ive  models  than  the  regular  solution  (and  quasichemical)  model  employed  by  Stringfellow. 

The  current  paper  attempts  to  assemble  a  data  base  for  III-V  compounds  which  can  be  employed 
to  compute  multicomponent  phase  diagrams  over  a  wide  range  of  temperatures (and  pressures)  both 
outside  and  within  the  quasi-binary  plane.  In  addition,  the  Cd-Hg,  Cd-Te  and  Hg-Te  binary 
systems  are  analyzed  and  combined  to  calculate  the  Cd-Te-Hg  ternary  system  over  a  range  of 
temperatures  (and  pressures)  in  order  to  compute  tie-line  compositions  between  quasi-binary 
CdTe-HgTe  alloys  and  the  Cd-Te-Hg  liquid.  The  former  alloys  are  of  current  interest  in 
synthesizing  controlled  band  gap  electro-optical  materials.  In  the  analyses  of  both  the 
III-V  and  1 1 -VI  (Cd-Te-Hg)  systems, the  methods  and  formalism  which  are  applied  are  those 
employed  in  dealing  with  metals  and  metallic  alloys  (12, 13),  silicides  (14),  and  oxynitride 
systems  (15). 

2.  Thermochemical  System  Employed  to  Characterize  Binary  III-V  Phase  Diagrams 

The  method  employed  for  describing  solution  and  compound  phases  is  the  same  as  that 
employed  to  describe  a  variety  of  metal,  metalloid,  and  oxide  systems  (13-15)  incorporating 

*  This  work  has  been  sponsored  by  the  Air  Force  Office  of  Scientific  Research,  Bolling  Air 
Force  Base,  Washington,  D.C.  under  Contract  F44620-80-0020  and  the  Electronic  Systems 
Division,  RADC-AFSC-Hanscom  Field,  Bedford,  Mass,  under  Contract  F19628-78-C-0060. 

On  leave  from  Department  of  Metallurgy,  University  of  Manchester-UMIST,  Manchester,  England. 
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|^some  sybolic  usage^which  facilitates  data  handling  as  indicated  below.  The  free  energy  of 
t--  the  liquid "phase.G  ,  in  the  binary  system  Al-In  is  given  by  equation  (1)  where  T  is  in 
[|  Kelvins,  R=8.314  J/g.at?K,  and  x  is  the  atomic  fraction  of  Indium:  1 

GL  =  ( 1 -x) °G^1 +x°G jn+RT(x£nx+ ( 1 -x)£n ( 1 - x) )  +x(l-x) [(l-x)LALIN+xLINAL] .  (1) 

i  Similarly*  the  free  energy  of  the  fee  phase  in  the  Al-In  is  defined  by  Equation  (2)  as 

XCC  r  CC* 

G  =  (1-x)°Ga1  +x°Gjn  +RT(x£nx+(l-x)  £n(l-x)>x(l-x) [ (l-x)AALIN+xAINAL]  (2) 

where  the  difference  in  free  energy  between  fee  and  liquid  aluminum, ALALIA,  and  the  difference 
in  free  energy  between  fee  and  liquid  indium,  ININLA,  are  listed  in  Table  1,  while  the 
^solution  parameters  LALIN,  LINAL,  AALIN,  andAINAL  are  given  in  Table  2.  TTie  liquid  para¬ 
meters  for  the  In-As  and  Ga-As  systems  were  taken  directly  from  Brebrick's  analysis  (11). 
Brebrick  also  provides  parameters  for  the  In-Sb  and  Ga-Sb  systems,  namely  LINSB=-3431-19.782T, 
LSBIN=-27698+19.581T,  LGASB=21652-32. 10ST  and  LSBGA=-30020+28. 259T  J./g.at.  While  these 
parameters  provide  a  good  description  of  the  In-Sb  and  Ga-Sb  phase  diagrams  and  the  excess 
free  energy  of  mixing  for  the  liquid  phases  in  these  systems, they  do  not  describe  the 
^enthalpy  of  mixing  for  the  liquid  phase  in  the  In-Sb  and  Ga-Sb  systems  as  measured  by  Ansara 
et  aL(7).  Thus  in  the  case  of  the  In-Sb  system,  the  enthalpy  of  mixing  defined  by  Brebrick's 
LINSB  and  LSBIN  parameters  is  given  by  Equation  (3) as 


ehl 


-x(l-x)  [(l-x)3431+27698x]  J./g.at. 


(3) 


(  where  x=Xg^.  By  contrast,  the  equation  for  EIIE  provided  by  Ansara  ot  ah  (7)  for  In-Sb,  which 
is  based  upon  experimental  measurements  of  the  heat  of  mixing,  is  given  by 

(4) 


EHL  -  -x(l-x)  ( 11689+1 9774 x-46661x2+24980x3)  J/g.at. 


Ansara  al. also  provide  the  following  equation  for  the  excess  entropy  of  mixing  of  the 
liquid  as 


ESL 


x(l-x) [4.602-7. 799x+12. 104x“~6. 782x  j 


J./g.at.  °K 


(5) 


Although  the  excess  free  energy  provided  by  Brebrick's  parameters  and  those  generated  by  Eqs. 

(4) and (5) in  the  temperature  range  of  interest  (near  700°K)aro  in  relatively  good  agreement, 
the  enthalpy  of  mixing  described  by  Equations (3)and (4) differ  by  800-1600  j/g.at.  Consequently. 
Equations  (4) and  (5) were  employed  to  derive  I.INSB=-15380-3.607T  and  LSBIN»-10293-2.393T  shown 
in  Table  2.  These  expressions  provide  a  closer  fit  to  Equations  (4 )and (5 )than  Brebrick's 
parameters  and  still  yield  a  satisfactory  agreement  in  the  calculated  and  observed  In-Sb 
phase  diagram.  A  similar  comparison  can  bo  made  for  the  Ga-Sb  system.  In  the  latter  case 
Brebrick's  analysis  yielded  LGASB=216S2-32. 108T  and  LSBGA=-30020+28.259T  J./g.at.  By  contrast, 
Ansara  et  al. (7)  suggest  that 


EHL  =  x(l-x) [531-21599x+32175x2-15788x3]  J/g.at. 
where  x^Xg^  and 

ESL  =  x(l-x)[4. 849-14. 4S6x+24.527x2-15.786x3]  J/g.at.°K 


(6) 


(7) 


The  present  parameters,  LGASB=4962-3. 209T  and  l.SBGA=-9715-0.456T  provide  a  better  description 
of  the  experimental  heat  of  mixing  for  Ga-Sb  alloys  embodied  (12)  by  Equation  (5)  than  do  the 
parameters  suggested  by  Brebrick.  The  parametric  descriptions  of  the  In-Al,  Ga-Al,and  In-Ga 
•  systems  shown  in  Tables  1  and  2  were  derived  from  the  earlier  description  provided  by 
Ansara  et_aL(12).  The  remaining  III  -  V  systems  were  characterized  as  shown  in  Tables  1-3 
based  on  the  available  thermochemical  and  phase  diagram  data  employing  the  standard  compila¬ 
tions  of  such  information  (16-20)  and  available  lattice  stability  data.  Table  4  compares 
the  calculated  and  observed  thermochemical  properties  of  III - V  compounds  based  on  this 
description  while  Figures  1-14  display  the  partial  phase  diagrams  computed  with  the  current 
^description.  The  heat  of  mixing  of  liquid  ,-Sb„  j.  at  1400 J K  is  given  as  -2477  J./g.at. 
based  on  experimental  data  (17)  and  -1126  J./g’.at.  off  the  basis  of  Table  2. 
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TABLE  1 


SUMMARY  OF  LATTICE  STABILITY  PARAMETERS 
(All  units  in  Joules  per  gram-atom  (mole  of  atoms),  T  in  Kelvins) 

L=liquid*,  A=fcc*,  T=tetragonal,  O=orthorhombic, 
R=rhombohedral,  X=red  phosphorus  structure,  W=white  phosphorus  structure 


ALALIA 

ALALLO 

Sagalo 

GAGALA 

GAGALT 

ININLT 

^NINLA 

ININLO 


4  *ALALLA 


10711  -  11.506T* 
0  -  18.451T 

5590  -  18. 45 IT 
-5021  -  8.368T 

-1674  -  8.368T 

3264  -  7.S94T 

2908  -  7.113T 

0  -  18.451T 


PPPPLX  «=  17991  -  21.506T 
PPPPWX  =  18171  -  19.539T 
PPPPLA  =  -18828  -  8.368T 
PPPPLR  =  4184  -  23.012T 

SBSBLR  =  19874  -  21.966T 
SBSBLA  *=  -18828  -  8.363T 

ASASLR  =  24874  -  23.012T 
ASASLA  =  -18828  -  8.368T 
ASASLX  =  21631  -  21.631T 


=  G^jqUld  -  G^C  =  10711  -  11.506  T  J./g.at. 


TABLE  2 

SUMMARY  OF  SOLUTION  PHASE  PARAMETERS* 

(All  units  in  Joules  per  gram-atom  (mole  of  atoms),  T  in  Kelvins) 

•*  *for  A£-In  system,  L= liquid,  A=fcc,  x=xJn 

GL=(l-x)Gj^  *  xGjn  ♦  RT(x£nx  ♦  (l-x)£n(l-x) )  ♦  x(l-x) [(l-x)LALIN  ♦  xLINAL] 

GfcC=(l-x)Gf»C  *  xG^C  ^  RT(x£nx  ♦  (l-x)£n(l-x) )  ♦  x(l-x)  [ (l-x)AALIN  *  xAINAL] 


LALAS  »  -16317  -  41.338T 
LASAL  =  -16317  -  41.338T 

LALPP  =  -38911  -  26.275T 
LPPAL  *  -38911  -  26.275T 

LINPP  *  -29288  *  5.188T 
LPPIN  -  -79496  ♦  5.188T 
RINPP  =  8368 
RPPIN  =  8368 

DGAPP  =  -90374  +  34.225T 
LPPGA  =  -58576  -  49.036T 
XGAPP  =  8368 
XPPGA  =  8368 

LGASB  =  4962  -  3.209T 
SSBGA  =  -9715  -  0.456T 
RSBGA  =  RGASB  =  8368 

LALGA  =  3280  -  2.021T 
LGAAL  =  2071  -  0.745T 
AALGA  =  -1046  +  5.439T 
IGAAL  =  -20962  +  26.778T 
OGAAL  =  OALGA  =  8368 


LALIN  =  25522  -  3.347T 
LINAL  =  20502  -  2.510T 
AALIN  =  60250  -  3.347T 
A INAL  =  55229  -  2.510T 

LALSB  =  8050  -  9.706T 
LSBAL  =  -17054  -  9.706T 

LSBAS  =  LASSB  =  0 
RSBAS  =  RASSB  =  8368 

LSBPP  =  LPPSB  =  -15899 
RSBPP  =  RPPSB  =  -20920-4. 184T 
XSBPP  =  XPPSB  =  -15899 

LASPP  =  LPPAS  =  -15899 
RASPP  =  RPPAS  =  -20920-4. 184T 
XASPP  =  XPPAS  =  -15899 

LGAAS  =  -1845  -  21.209T 
LASGA  =  -41551  +  7.192T 
RGAAS  =  RASGA  =  8368 

LINAS  =  -243  -  29.522T 
LASIN  =  -48124  +  17.472T 
TASIN  =  TINAS  =  8368 


TABLE  2  (CONTINUED) 


SUMMARY  OF  SOLUTION  PHASE  PARAMETERS 
(All  units  in  Joules  per  gram-atom  (mole  of  atoms)  T  in  Kelvins) 


J.INSB  =  -15380  -  3.607T 
I.SBIN  =  -10293  -  2.393T 
I  TSBIN  »  TINSB  =  8368 
RSBIN  *  RINSB  =  8368 


LGAIN  =  4435  +  2.268T 
LINGA  =  4435  +  0.954T 
TGAIN  =  -3431  +  21.907T 
TINGA  =  -3376  +  24.393T 
OGAIN  =  OINGA  *  12552 


TABLE  3 


GS_  q  j 

i 


Compound 


COMPOUND  PARAMETERS  FOR  III-V  ZINC-BLENDE  COMPOUNDS 
+  0.5  Gjc  +  0.25  (0.5  LIIJJ  +  0.5  LJJII  -  C]  Joules/g.at. 


:  "  Ga0.5As0.5 


Compound  Parameter  (C)  Compound 
(Jouies7g.at.  bK) 


259,525  -  83.496T 


Compound  Parameter  (C) 
(Joules/g.at.  **10 


A£0.5Po.5 

364,443  - 

60.  517T 

Ga0.5Sb0.5 

179,343  - 

77.571T 

1  A£0.5AS0.5 

305,348  - 

78.994T 

In0.5P0.5 

188,749  - 

64.099T 

^o.s^o.s 

174,992  - 

55.815T 

In0.5As0.5 

187,217  - 

73.806T 

Ga0.5P0.5 

280,328  - 

75.647T 

In0.5Sb0.5 

125,227  - 

63. 011T 

COMPARISON 


* 

! 

r 

;  Compound 


TABLE  4 

OF  CALCULATED  AND  EXPERIMENTAL  THERMOCHEMICAL  PROPERTIES  OF  III-V  COMPOUNDS 

(J./mol . ,  J./mol . °K) 

AHf[298°K]  ASf[298°K]  AS  [fusion] 


calculated 

exptl.  calculated 

exptl. (16) 

calculated 

exptl. 

\lP 

-164,849 

-164,431+2929(16) 

-3.97 

-3. 85i5. 0 

61.66 

61.50*1.6(16) 

MAs 

-117,127 

-122,591*5021(16) 

-4.18 

-3.76+4.8 

70.91 

70.42(43) 

MSb 

-51,045 

-50,208*4 18(16; 

-9.46 

-8. 86+0. 8 

59.32 

61.67(43) 

GaP 

-129,971 

-12^173+8368(16) 

-10.92 

-11. SI*  1.9 

66.10 

66.53(43) 

GaAs 

-f‘6,291 

-88,575(11) 

-15.51 

-12.47+1.3 

70.015 

69.63(11) 

GaSb 

-41,547 

-41,840*1255(16) 

-14.19 

-9.20+1.9 

67.061 

67.31(11) 

InP 

-84,391 

-75,312*8368(16) 

-21.02 

-20.88+1.9 

59.07 

59.08(43) 

InAs 

-61,647 

-61^23(11) 

-18.74 

-18.83+2.5 

63.92 

63.36*6.89(16) 

InSb 

-29,974 

-30,962(11) 

-15,93 

-15.69+1.1 

58.53 

59.82(11) 
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L  (liquid) 


Figure  1.  Calculated  (la- In  Phase  Diagram 


Figure  2. 


Calculated  In  -  Sl>  Phase  Diagram 
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Figure  3.  Calculated  Al-Ga  Phase  Diagram 
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Figure  4.  Calculated  In  -  As  Phase  Diagram 
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The  free  energy  of  the  ternary  liquid  phase  in  the  P-In-As  system  is  defined  on  the 
basis  of  Kohler's  equation  as 

GL  =  (l-x-y)0Gp+x°(ij;I'+y0G^s+RT((l-x-y)tn(l-x-y)+x!nx+yiny) 

♦  O-x-ylxCl-yf^O-x-yJlPPIN+xLINPP)  +  (l-x-y)y(l-xf 1  [(l-x-y)LPPAS+yUSPP) 

+  xy (x+y) " * [xLINAS+yLASIN]  *  xy(l-x-y)  TR.NL  J./g.at.  (8) 

where  x  is  the  atom  fraction  of  In,  y  is  the  atom  fraction  of  As  and  TRNL  is  the  ternary  inter¬ 
action  parameter.  Values  of  TRN1.  for  five  ternary  systems  which  have  been  investigated  here 
are  listed  in  Table  5.  The  specific  value  for  the  P-In-As  system  is  12552  J/g.nt.  This  is  a 
rather  small  value  since  the  maximum  value  of  the  ternary  term  occurs  at  x=y* (l-x-y)*l/3.  Thus 
the  maximum  contsibution  of  this  term  is  only  465  J/g.nt.  If  attention  were  restricted  to  the 
quasi-binary  join  (i.e.,  where 'X=0.5)  between  In-  ,1’  -  and  Inn  -As-  -  the  maximum  contribution 
is  390  J/g,at.  when  x»0.5,  y»0.2S  and  (l-x-y>0.2SV  1  ’ 

The  free  energy  of  the  In-  ,!’  ,-In-  rAs„  ,  zinc -blende  phase  (S)  is  defined  by  Equa¬ 
tion  (9)  ns  follows  0-5  1,5  '•J 

GS  =  (0.5-y)°GpCC  ♦0.5oc[te(C*yoG^c  •*  (l-2y)  (0.25)  fO. f»  I.PPlN»O.SI.!NPP-Cj) 

♦  2y (0.25)  (0.5I.INAS*0. 5EASIN-C ,)  *  CAB  y(l-2y) 

»RT(ytny* (0.5-y) tn(0.S-y)-0.5!nO.S)  J./g.at.  (9) 

where  Cj.  the  compound  parameter  for  In^  ,.P()  ^  =  188,749-64 .099T  J/g.at.  and  C^,  the  compound 
parameter  for  Infl  ^As^  187,217  -  73.806T  J./g.at., as  given  in  Table  3.  The  interaction  para¬ 
meter  for  mixing  ln^  jP^  ^  and  In^  -As^  CAB,  is  equal  to  zero  in  this  case,  as  shown  in 

Table  5.  Thus  this  compound  is  defined  as  an  ideal  solution  of  its  components.  Since  the  free 
energies  of  the  liquid  and  zinc-blende  phases  can  be  defined  explicitly,  the  partial  ternary 
phase  diagrams  covering  tho  I.-S  equilibria  can  be  computed  as  a  function  af  temperature.  This 
procedure  has  been  carried  out  for  the  ternary  systems  listed  in  Table  5.  'Hie  results,  dis¬ 
played  as  isothermal  sections  or  as  quasi -binary  joins. are  compared  with  experimental  findings 
and  Stringfellow' s  calculations  (1)  in  figures  15-25.  A  calculated  Isothermal  section  due  to 
Ansara  et  aj,  (7)  for  the  Sb-Ga-ln  case  at  825“)’ is  included  for  comparison  in  Figure  21.  Tho 
latter  results  are  in  good  agreement  with  the  current  findings  shown  in  Figure  2  3,  Although 
it  is  not  possible  to  make  an  exact  comparison  of  the  interaction  parameters  employed  by 
Stringfellow  in  his  quasi-binary  calculat  ions  (1,6)  and  iho.,e  employed  here,  an  approximate  com¬ 
parison  can  be  made  in  the  following  way.  In  Stringfcllow's  study. the  quasi-binary  system 
InP-lnAs  was  treated  as  if  the  liquid  and  solid  phases  were  both  regular  solutions.  The 
interaction  parameter  for  the  liquid  was  taken  to  be  5238  J^mol.  Since  the  heat  of  mixing  is 
the  product  of  the  mol.  fraction  and  the  interaction  parameter,  the  maximum  value  of  the  heat  of 
nixing  is  (0. 25) (5238)  =  1309  J/mol*65$  J/g.at,  for  liquid  0.51n0  j-O.Sliig  ..As^  j.  This 

value  is  listed  in  Table  6  along  with  comparable  parameters  for  tho  other  systems  of  current 
interest  cited  by  Stringfellow.  Tho  latter  results  are  compared  with  those  derived  from  the 
current  study  on  the  basis  of  Equations  (8)  and  (9)  and  Table  5.  As  indicated  above, the  maxi¬ 
mum  contribution  along  the  quasi-binary  join  is  equal  to  TRNI./32  and  CAB/8.  Table  6  lists 
these  values  for  comparison  with  Stringfcllow's  results  in  Table  6.  While  there  is  no  obvious 
correlation  between  the  current  results  and  Stringfcllow's  values, all.  of  the  values  listed  are 
small!  The  analysis  of  the  Ga-Sb-In  system  carried  out  by  Ansara  et  al.  (7)  utilized  a  ternary 
heat  of  mixing  equal  to 

EIIL  3  xInxGa*Sb  (23234x,  -24066xr  *582x_.  J -  (10) 

1  In  Ga  Sb’ 

The  present  analysis  sets 

F.  L  *  -x.  x,.  X-.  8368  .  (11) 

it  In  ua  Sb 

At  x.  =x_  *x-,=0. 333, Equation  (10)  yields-3.1  J/g.at. .while  Equation  (11)  yields  -309  J./g.at. 
At  coSpositions  near  the  edges  of  the  ternary  (0. 1 ,0.4 ,9.5)  Eq.  (11)  yields  -167  J/g.at. .while 
Eq. (10)  yields  results  between  +138  and  -188  J.g./at.  A11  of  these  contributions  are  small. 
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TABLE  5 


COMPILATION  OF  TERNARY  PARAMETERS  FOR  III-V  LIQUID  AND  ZINC-BLENDE 

COMPOUND  PHASES 

(All  Units  in  Houles  per  gram -atom  (mole  of  atoms)  T  in  Kelvins) 


System 

TRNL 

(Joules/g.at.) 

Zinc-blendc  Phase 

CAB 

(Joules/g 

P-In-As 

12552 

In0 .5(P,As)o.5 

0 

As-In-Sb 

0 

In0.5^As,Sb^0.5 

9623 

P-Ga-As 

-31380 

Ga0.5(‘P,As^D.5 

■  2092 

Ga-Sb-In 

-8368 

(Ga,In)o5Sb0.5 

4602 

Al-Ga-Sb 

-8368 

(A1’Ga)0.SSb0.5 

-4184 

TABLE  6 

COMPARISON  Of  THE  MAXIMUM  QUASI-BINARY  HEAT  OF  MIXING  FOR  THE  LIQUID  AND  SOLID  COMPOSITIONS 
FROM  STRINGFELLOW  (1,6)  AND  THE  CURRENT  ANALYSTS  (Joules/g.at.) 


System 


LIQUID 

Stringfel low(l)  Current 


SOLID 

Stringfel low  Current 


In0.SP0.5-In0.S  AS0.5 
,n0.5As0.S'In0.5Sb0.5 
(,a0.5l>0.S’(,a0.SAs0.5 


Gn0.SSb0.5“In0.5Sb0.S 

A*0.5Sb0.5'C,O0.5Sb0.S 


Thus  the  current  analysis  suggests  a  simple  method  for  calculating  isothermal 
■Sections  for  1 1 1 - 1 V  ternary  systems.  In  particular,  if  a  given  section  is  required,  the  data 
base  contained  in  Ta hies  1-3  can  be  employed  to  calculate  the  section  of  interest  with  TRNL= 
QAIS=0.  The  process  can  then  be  repeated  employing  a  range  of  TRNL  and  CAB  similar  to  those 
Shown  in  Table  5  to  predict  the  spread  of  expected  results.  If  a  few  experimental  data  are 
mailable  at  a  given  temperature,  they  can  be  readily  employed  to  fix  TRNL  and  CAB,  which  in 
-jnirn  can  be  employed  to  predict  the  behavior  at  other  temperatures.  It  is  hoped  that  future 
,  studies  will  permit  improvement  of  the  current  data  base. 

3.  Thermochemical  System  Employed  to  Characterize  Binary  II-VI  Phase  Diagrams 

The  1 1 - IV  compounds  of  Cd„  ^Te^  r  and  llgQ  5TeQ  ,  form  a  series  of  cubic  zinc-blende 
solid  solutions  which  offer  the  attractive  property'of  a’variable  band  gap  that  is  a  function 
,  composition.  The  (Cd,Hg)0  5TeQ  ^solid  solution  displays  band  gaps  which  run  in  a  nearly 

I  inear  fashion  between  the  wide  gap  semiconductor  CdQ  jToq  (E^  =  1.6  eV)  and  the  semi- 
metallic  compound  HgQ  ^TeQ  s>which  can  be  considered  to  be'a  semiconductor  with  a  negative 
band  gap  (37).  Alloys  of  Cd0  lHgQ  4TeQ  $  and  CdQ  2Ug0< jTeQ  exhibiting  band  gaps  ranging 

■from  0.1  to  0.5  eV,  respectively,  are  of  particular  interest  for  detection  of  infrared  radia- 
-j-ion  in  the  2  to  20  micron  range.  The  growth  of  such  alloy  compositions  from  the  melt  presents 
*  Some  difficulties  due  to  the  "width"  of  the  two  phase  liquid  plus  solid  field  along  the 


1323  I-  O  ref.  (20) 
□  ref.  (21) 
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Figure  15.  Calculated  and  Observed  Equilibria  Along  the 
ln0.5As0.5  '  ,n().5i’().5  -binary  Join 
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Figure  16.  Calculated  and  Observed  Equilibria  Along  the 
In0.5Sh0  5  *  Ino  SAs0  5  Q»asi-binary  Join 
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Calculated  Partial  Isothermal  Sections  in  the  Phosphorus 


Calculated  Partial  Isothermal  Sections  in  the  Arsenic- Indium-Antimony  System 


sothcrmal  Sections  in  the  Gal  1  ium-Ant imony- Indium  System 


CdO.STe0.5'tlg0.5Te0.S  ^uasi'binary  join.  Liquid  phase  epitaxial  (LPE)  growth  methods  have 
been  employed  to  circumvent  this  problem  (37),  but  controlled  growth  requires  a  knowledge  of 
the  tie-lines  between  the  alloy  compound  and  the  liquid.  Such  information  can  be  generated 
experimentally  with  considerable  effort.  Alternatively,  it  is  possible  to  use  the  CALPHAD 
method,  which  has  been  applied  to  a  very  wide  range  of  systems,  for  calculating  the  ternary 
Cd-Te-Hg  system  from  a  knowledge  of  the  component  binary  systems  in  order  to  provide  the 
required  information  on  tie-line  composition.  The  description  provided  below  provides  an 
account  of  such  calculations. 

4.  Description  of  the  Vapor,  liquid  and  Zinc-blende  Phases 

Table  7  summarizes  the  current  descriptions  of  the  liquid(L),  vapor(V)  and  zinc- 
blende(S),  phases  in  the  cadmium-tellurium-mcrcury  system  as  a  function  of  temperature  (T,°K), 
pressure  (P,  atmospheres),  atomic  fraction  tellurium(x) ,  and  atomic  fraction  of  mercury (y) . 
These  equations,  along  with  the  associated  lattice  stability  values  relating  the  free  energy 
of  the  vapor,  liquid,  fee,  and  stable  forms  of  cadmium,  tellurium  and  mercury, have  been  assem¬ 
bled  by  employing  oxisting  (38,14-17)  thermochemical  and  phase  diagram  data  for  the  pure 
elements  and  binary  systems.  Figures  26-31  show  the  partial  binary  phase  diagrams  computed 
for  the  Te-Cd,  Hg-Cd  and  Te-Hg  systems  at  P=1  atmosphere  and  P«16  atmospheres.  Thcso  phase 
diagrams  follow  directly  from  the  ternary  equations  listed  in  Table  7  when  they  are  reduced 
to  binary  systoms  for  each  of  the  edgo  binary  systems  in  question.  The  small  ternary  terms 
assigned  to  the  liquid  phase  (-33472xy (1-x-y) )  and  the  zinc-blonde  phase  (-4180y  (1 -2y) ) 
Joulos/g.at.  were  chosen  in  conformity  with  tho  experimental  rosults  along  the  CdQ  ,.To0  j- 

l*Ko  s^eo  5  quasi-binary  join  shown  in  Figures  32  and  33.  Sinco  the  maximum  valuo  of  these 

terms  is  attained  at  x=0.S  and  y*0.25,  it  is  apparent  that  the  ternary  correction  term  is  Bt 
most  about  -1000  Joules/g.at.  for  the  liquid  phase  and  -500  Joules/g.at,  for  the  zinc-blendo 
phase  at  the  x=0,5,  y»0.2S  composition.  Table  8  compares  calculated  and  observed  thermocheai* 
cal  values  for  Cd0>5Te05  and  llgo.sTco.S  at  298  K. 

Figures  34-36  display  calculated  isothermal  sections  in  tho  Cd-To-Hg  systems  between 
1  and  74  atmospheres  at  temporaturcs  ranging  from  773°K  to  1213<>K.  Those  calculated  iso¬ 
thermal  soctions  show  the  tie-lines  between  tho  zinc-blende  and  liquid  phases  in  addition  to 
the  three  phase  vapor /liquid/zinc-blcnde  fields  which  are  especially  important  in  crystal 
growth  operations.  The  calculations  of  tio-lines  can  be  compared- directly  with  the  recent 
experimental  rosults  of  Mroczkowski  and  Vydynnath,  who  determined  tie-line  compositions  by 
means  of  "closed  tubo  tipping  experiments".  In  those  experiments,  crystals  of  the  zinc-blende 
phase  of  fixed  composition  are  equilibrated  with  the  liquid  phase  in  a  graphite  boat.  After 
equilibration,  the  temperature  was  decreased  by  5  to  7°C  and  the  liquid  was  decanted.  Subse¬ 
quent  eloctron  microprobo  analysis  of  the  concentration  profiles  along  the  length  of  the 
crystal  (from  the  zinc-blendo  single  crystal  phase  into  the  polycrystallino  phase,  which  forms 
when  the  contact  liquid  freozes)  permits  determination  of  tho  tie- line  compositions  (37).  At 
818°K,  the  following  comparison  can  be  made  between  the  calculated  and  observed  tie  lines. 

Zinc-blende  Liquid 


Atomic  Percent 

!!& 

Cd 

!!& 

Cd 

i 

Observed  (37) 

32.5 

17. S 

20.3 

2.3 

Calculated 

32.5 

17.5 

19.5 

2.0 

i 

This  excellent  comparison  provides  an  independent  check  of  the  description  of  the 
liquid  and  zinc-blende  free  energies  shown  in  Table  7.  These  equations  can  be  used  to  com¬ 
pute  the  liquid/vapor/zinc-blende  equilibria  over  a  wide  range  of  temperatures  and  pressures, 
which,  in  turn,  can  be  used  to  define  the  conditions  for  stable  conventional  and  LPE  growth  of 
Te-rieh  and  Tc-poor  (Cd.Hg)^*^  eorap05itions. 

5.  Conclusions 

The  examples  of  calculated  III -V  and  II-Vl  multi -component  systems  presented  here 
illustrates  the  scope  and  utility  of  the  method  presented.  The  range  of  systems  to  which  such 
computations  can  be  applied  will  expand  as  additional  binary  systems  are  added  to  the  data  base. 
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TABLE  7 


SUMMARY  OF  FREE  ENERGY  EQUATIONS  DESCRIBING  TOE  VAPOR 
LIQUID, AND  ZINC-BLENDE  PHASES  IN  TOE  CADMIUM-TELLURIUM-MERCURY  SYSTEM 
(All  units  in  Joules  per  gram-atom  (mole  of  atoms),  °K) 

L= liquid  phase,  V=vapor  phase,  S=zinc-blende  phase,  A=fcc 
R=rhombohedral  (stable  Hg) ,  T=trigonal  (stable  Te) ,  E=hcp  (stable  Cd) 
x=atom  fraction  Te,  y=atom  fraction  Hg,  P=pressure  (atmospheres) 

R=8.314  J/g.atfK 


Lattice  Stability  Values 


G*J  -  G^  =  2292-9. 832T;  T=234°K 

Clfi  -  C»5  =  2297-1.464T 

Hg  Hg 

g}'  -  gJ:  n  59174-93, 303T+RT£nP 

Hg  Hg 

G,p  -  Gl!  =  1 7489-24. 184T;  T=  723°K 
le  Le 

G.J  -  cl  =  17489-15. 816T 

!c  Te 

'■fil  -  G.p  =  58409-45 . 982T+0 , 5RT£nP 
le^  ie 


GCd  "  4  =  6192-10-418T;  t=594K 

gJ.  -  g!.  =  1924-2. 050T 
Cd  Cd 

.  -  cj;.  =  102, 520-98. 659T*RT£nP 
Cd  Cd 


Ternary  Phases 

GL  =  (l-x-y)G^  *  xG^e  ♦  ycjj  +  RT  [ (l-x-y)£n(l-x-y)+x£nx+y£ny] 

-xy(x+y)  1  [( 30962*5. 02T)x*( 35 14*1 6. 74T)y)  ♦  y(l-x-y) [-10878*4 . 184T] 

♦x(l-x-y) (1-y)’1 [ (1-x-y) (112968-159T)  -x(75312*20.92T)) -  33472  xy(l-x-y) 

V  V  V  V 

av  .  (i-x.y)G*d  .  o.sx  ♦  >• 

♦  RT((l-x-y)£n(l-x-y) +x£nx*y£ny) 

GS  =  (0.5-y)  G^j  ♦  0.5  *  y  G^  *  (1 -2y) [- 61 19 1 ♦ 13.075T] +2y[-26798*7. 11T] 

+  RT[y£ny  +  (0.S-y)£n(0.5-y)-0.5£n0.5)  -4180y(l-2y) 

0  <  y  <  0.5 


TABLE  8 

COMPARISON  OF  CALCULATED  AND  OBSERVED  (16) FREE  ENERGY  OF  FORMATION  AT  298°K 


Compound 


Cd0.Sre0.5 


Hg0.5Te0.5 


AH  [2981  (J./g.at.) 

-50900  ±  400  (obs.) 

-51500  (calc.) 


-15900  t  2100  (obs) 
•8""  '  * 


AS  [2981  (J./g.at.°K) 

-4.18  +  1.25  (obs.) 
-4.14  (calc.) 

-6.28  ±  1.26  (obs.) 


Molt  Fraction  x 

•  -  Steininger  (39) 

0  Uouidut  Data  -  Blair  iryj  Hewnham  (40) 

V  Solxlut  Oata  -  Harrnan  and  Strauu  (4 
□  Solidvu  Data  -  Schmit  and  Spttrichntidtr  (42) 

Figure  32.  Liquidus  and  Solidus  Curves  for  Pseudo 
binary  Hg^_x^CdxTo  Compositions  (39) 
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Figure  33.  Composite  P-T-x  Diagram  for  the  Pseudo¬ 
binary  HgTe-CdTe  System  (39) 


Calculated  Isothermal  Sections 
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Figure  33.  Calculated  Isothermal  Sections  in  the  Cd-Te-lte  Svitpm 
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Figure  36.  Calculated  Isothermal  Sections  in  the  Cd-Te-Hg  System  at  1213°K 
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CALCULATION  OF  ORDERING  TEMPERATURES  FOR  THE 
QUASI-BINARY  SECTION  Fe  Al  -  Fe  Si 


1.  Introduction 

Where  limited  experimental  data  is  available,  the  method  adopted  by  Inden  (1-7) 
is  attractive  because  it  requires  only  two  ordering  parameters  for  each  of  the  three  component 
binary  systems,  provided  the  temperatures  of  interest  are  above  any  magnetic  transition 
temperatures.  Those  ordering  parameters  are  the  first  and  second  nearest  neighbor  (chemical) 

(1)  p) 

interaction  parameters  denoted  Wv  and  W  .  These  can  be  related  to  heats  of  formation  as 
well  as  critical  ordering  temperatures,  as  will  be  shown  in  the  following  section.  The  Fe-Si 
binary  system  (Fig.  1)  has  already  been  characterized  in  some  detail  (3,4)  and  the  previously 

determined  values  of  K ^ 1  ^ :mci  W^Fe-Si  ^avo  ,)0Cn  useci  unaltered  in  the  present  work 

Table  I).  Their  validity  has  been  confirmed  by  the  successful  prediction  of  ordering 
equilibria  in  ternary  Fe-Si-Co  alloys  (8).  The  Fe-Al  system  is  inherently  more  complicated 
19,10)  (Fig.  2)  and  several  additional  effects  have  been  recently  added  by  KUster  et  al  (11). 
Nonetheless  the  major  features  have  been  reproduced  by  t he  simplified  Inden  model (12).  Only 

a  minor  change  has  been  made  in  the  value  of  1  ^ j.'o- Al  usetl  l>rcv*ously  by  Inden  (6)  in  the 

light  of  the  available  experimental  data  for  Fe-Si-Al  alloys  (Table  II).  The  derivation  of 

(11  P) 

suitable  values  of  suitable  Wv  si_A1  and  K  Si_,u  values  poses  a  problem  in  so  far  that  all 

experimentally  accessible  properties  for  these  alloys  refer  to  either  FCC  Aluminum  or  diamond 

cubic  Silicon.  One  has  therefore  to  extract  the  relevant  values  of  h'O)  and 
,  ^  a  S  i  -  A 1 

K  A1  from  whatever  data  is  available  for  BCG  ternarv  alloys, 

Si-Al  •  1 

2 .  Experimental  Information  on  Fe-Al -Si  Alloys 

Lihl,  Burger,  Sturm  and  libel  (14)  have  given  a  general  indication  of  the  location 
of  the  ordering  region  in  Fe-Al-Si  alloys  (Fig.  3).  However  there  is  considerable  divergence 
between  this  estimate  and  the  investigation  of  Katsnel'son  and  Polishchuk  (13).  More  weight 
has  been  given  to  the  latter  work  as  it  utilized  relative  intensities  of  superlattice  lines 
(Fig.  4)  whereas  the  work  by  Ebel  et  al  merely  interpreted  changes  in  d- spacing  as  ordering 
effects.  An  NMR  study  of  Fe.,Si  .Al  alloys  (15)  suggests  that  Aluminum  and  Silicon 

l  *  **  K  J  X 

substitute  interchangeably  on  one  specific  sublattice,  which  confirms  that  the  interaction 
energy  between  Silicon  and  Aluminum  is  relatively  low. 


/  A 


\  xtk//' 

y  L'\!  A// 

rm 


ri  N  *  >'  / 

i  '  i  <  V 

I  X  I  N 

.  J* 

|  ;  /  '' 


4‘  !  * 


1  tl  \\^  ^ J  j  ^O;..  ||  } 

../  «to  i^waarascansJ*'  < 

<%C  >W*<  X.* 

V  7/  ••~,^> ^  X  . 


. . ^  X  ,  4 

(C>  '""^'"•••■'.Nl  '-SV 


Fig.  3  Suggested  Ternary  Phase  Boundaries 
for  Fe-Al-Si  alloys  from  F.  l.ihl  et  al  (14) 
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l'ig.  4  Critical  Temperatures  suggested  by  the  temperature 
dependence  of  superstructure  lines  from  Katsnel'son  (13) 
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3.  Ordering  Model 


(a)  Definition  of  order  parameters 


The  main  feature  of  the  BWG  model  is  the  description  of  atomic  distributions  inside 
specifically  defined  sublattices.  Thechoice  follows  from  the  hind  of  interaction  between 
atoms.  Any  BCC  lattice  may  be  divided  into  four  FCC  sublattices,  I,  II,  III,  and  IV,  with 

the  atomic  configurations  described  by  occupation  probabilities  ' p . ^ '  of  sublattice  site  ’ L* 
by  the  component  ’ i ' .  (Fig.  5)  1 

Sites  of  sub-lattice  I  — o 

Sites  of  sub-lattice  II—®  /  / 

'  Sitesof  sub  lattice  III—  A  /  /  / 

Sites  of  su b  lattice  I V— A  r/_ _  /  y 


/U~ j  ~vy 

~mrn 
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Fig.  5  Relative  position  of  sub  1: it t ices  in  the  DO „  structure. 

The  atomic  configuration  in  a  binary  alloy  A(,Bj_(.  can  then  be  dascuVd  with  the  aid  of  the 

following  three  parameters,  which  are  a  combination  of  the  independent  occupation  probabilitie 

L. 

X  =  1  (PA‘  *  p”  -  PAIn  -  PArV)  (la) 

y  =  l  (pAni  -  pAIV)  ot) 

z  -  1  Cpa1  -  PAlr)  (1c) 

) 

The  greatest  value  of  these  parameters  corresponds  to  the  highest  degree  of  order,  and  the 
concentrations  of  A  and  B,  (C^  and  C^J  set  limits  on  the  values  of  x,  y  and  z.  For  the 
binary  case  these  are  given  by: 
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(a) 

0S< 

x  j  CB 

for  0  <  Cg  <  0.5 

(2a) 

0>) 

0S< 

y  ^  min 

(CA  -  x,  CB  *  x) 

(2b) 

(c) 

0s< 

z  <  min 

(CA  *  x>  CB  " 

(2c) 

Case  1 


I'he  random  BCC  lattice  (A2)  has  an  equal  distribution  of  the  atoms  on  the  four  sublattices 
and  corresponds  tox=y=z=0. 

Case  2 

The  state  x/0,  y=z=0  indicates  'A'  preferentially  occupying  sublattices  I  and  ! I  (x  >  0) , 
or  III  and  IV  (x  0).  It  follows  that  the  number  of  AB  neighbors  in  nearest  neighbor  (n.n) 
positions  is  increased,  i.e.  this  corresponds  to  the  B2  structure. 

Case  3 


It  follows  that  y  *  0,  describes  a  surplus  of  A  atoms  in  next  nearest  neighbor  (n.n.n) 
positions,  which  corresponds  to  t he  IU)„  structure.  As  far  as  t he  present  problem  is  con¬ 
cerned  we  need  not  consider  z.  (See  I'efs.  1,2) 


( 1) j  free  lincrgy  repressions  and  Interaction  Parameters 


The  most  stable  atomic  configuration  then  follows  from  the  minimum  value  of  the 


.BCC 


.•onfigurational  free  energy/mole  AC  ‘  which  is  given  by: 


rBCC 

Mi 


BINARY 


=  [J°  -  NC^Cjj  (4W11)  ♦  3Wi?))  - 

+  31V C 2 )  (y-  f  2* ) 


N 


(8  w(!)  -  (M2))  x2 


NkT  L* 


~T~  i. 


L,  l  s  Llr  x  L 
PA  In  P.\  r  Pj)  ln  Pb 


(3) 


IJU  is  t tie  internal  energ>  of  puic  components;  the  second  term  is  the  energy  of  mixing;  tiie 


thud  tcim  gives  the  energy  contribution  due  to  ordering,  and  the  fourth  term  is  the  entropy 
contribution  of  the  configuration  which  is  also  expressed  in  terms  of  x,  y,  :.  N  is  Avagadro’s 

number,  and  k  the  Boltzmann  constant.  Iv  and  >\  '  are  defined  below  and  lor  eqn  (3)  to  be 


•3 


-24  -24 

valid,  these  must  be  expressed  in  k  units  where  1  k  unit  =  3.3  x  10  cal  or  13. S  x  10  J. 


Ihe  \alues  of  and  are  lelated  tu  tin.  nearest  and  next  nearest  neiglibor  bond  energies, 


1 1 


^  as  follows: 


K(0  -  -  2VW<‘>  *  VM(»  *  VBB0) 


w(2)  "  -  -vab(2)  ’  vo\(2>  *  %(2) 


(4a) 

(4b) 


According  to  this  definition,  positive  values  for  and  correspond  to  greater 


affinity  between  unlike  atoms.  All  interchange  energies  are  presumed  to  be  independent  of 
tempera tin e  and  of  environment  of  other  atom  sites  where  the  interchange  process  is  considered. 
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(11  (21  BCC 

If  the  values  of  W  and  W  are  known,  AG  '  can  be  calculated  with  respect  to  the  ordering 
parameters  x,  y,  2  using  the  necessary  equilibrium  conditions  including  particularly: 


3G  3G  3G 

3X  Dy  3; 

and  other  subsidiaiy  conditions. 


■  (t.  c) 


( c )  Expressions  for  critical  ordering  temperatures  in  binary  alloys 

The  critical  temperatures  T  (A2/H2),  T  (!52/DO_),  define  the  regions  in  which 

minimum  free  energy  is  obtained  by  the  A2,  112  and  DO,  configurations  respectively.  For  the 
binary  case  these  are  given  by: 

kT  =  kfV2/R2  =  (SiVU)  -  «*(?))  .  c  (• 

x  (3WG)  Bh’G  BMC  ■'  1 


kT  „  kT82/no3  =  ()lv(2)  .  (c  _  x)  cc  +  x) 


(BWG) 


wliere  x  is  1 110  appropriate  value  of  the  ordering  mot  ficieiit  at  the  temperature  T  . 

Since  lr  '  are  energy  parameters,  the>  are  most  oasil>  determined  from  energy  measurements 
like  enthalpies  of  mixing  of  random  alloys,  or  enthalpies  of  formation  All  (C,T)  of  alloys 
with  given  atomic  configurations.  These  entities  can  immediately  bo  expressed  by  means  of 

IV ^  for  a  binary  alloy  Aj  (,U,.  with  role  mice  to  the  nine  components  in  the  same  crystal 

structure  as  the  alloy  formed,  lor  !U‘G  alloys: 


A  2 

All 

»  -  N  C(1  -  C)  j  4WU)  a  5wU) 

....  (7) 

ah»2 

»  -  N  J^Cl-lwCO  a  3W^):-  dC2  h’(2)  j 

for  C  <0.5 

ah  1X33 

=  -  N  C  (4W(0  +  5W^2)) 

for  0  :  C  <  0.2S 

OR 

-  -  N  l.SW^)  a  C(4WC0  -  3W(2))  j 

for  0.25  <  C  <  0.5 

AHR32 

•—  «-J 

=  -  N  f  (4A'(  ‘-I  a  :>!*('))  -  4C“  :vlO  1 

for  C  <  0.5 

values  obtained  from  an  analysis  of  critical  temperatures  (e.g.  via  equations  oa,  6b), 
HlSu  .  ^ 

can  be  related  to  t he  K  values  ii*.  thermochcmieul  equations  (e.g.  equ.  71  by  a  scale  factor 
x  (listed  in  Table  1).  See  References  (11,2).  (Hie  differences  between  1\’W  and 

*1*)  are  essentially  due  to  the  omission  of  short  range  order  in  the  BKG  model.) 

BIVG 
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*  4 .  Expressions  for  the  Free  Energy  and  critical  ordering  temperature  in 
ternary  alloys  as  used  in  the  present  work 

In  the  ternarv  case  there  is  a  slight  difference  in  the  definition  of  the  order 
parameters  (5).  If  the  total  number  of  atoms  is  N,  distributed  over  4Nq  lattice  sites,  the 

occupation  probabilities  p^  (i  =  A,B^C,;  L  =  I,  II,  III,  IV)  are  combined  to  give  the  order 

parameters  as  follows: 

(al  Xj  -  i  (Pj1  *  PjI!  '  Pj"1  -  PjIV)  f*a> 

(b)  •  !  (P^11  -  P-jIV)  (8b> 

(=)  z}  -  1  (pd  -  Pj'")  l8c> 

with  j  -  A,  B 

A  similar  condition  to  Equation  (5)  occurs  except,  now  there  arc  two  x  parameters,  two  y 
parameters  and  two  z  parameters,  corresponding  to  the  degree  of  order  of  two  solutes  and  not 
just  one.  For  the  purpose  of  determining  and  T  temperatures  it  is  only  necessary  to 

consider  and  xR.  The  limiting  values  of  x .  ,xR  in  the  ternary  case  are: 

fa)  The  maximum  value  of  xA  is  the  smaller  of  a  or  (1-a).  (9a) 

(b)  The  maximum  value  of  xR  ir  the  smaller  of  b  or  fc  •  xA)  (9b) 

In  the  specific  case  of  Fc,Al ^ ^  ^  S  i  ^  alloys,  where  a  1  Fc,  b  -  A1  and  c  =  Si ,  this  leads 

to  a  situation  where  (a)  is  always  0.75  and  therefore  x.,  .  =  0.25.  Likewise  since 

A  (.  ma  x ) 

c  =  (1-a-b),  x.  ,  ,  =  -b.  The  configurational  free  energy/mole  All  of  the  three-component 

iHmaxj 

solid  solution  is  then  given  by: 


ac(bcc)  -  u0  -  N  (4  CCAcBwCO^  ♦  W(')AC  *  cB.;cw(‘'k  ) 

+  3(caCbw(2).\B  +  CACCW(2)AC  +  CBCCly(2)BC  ^  } 

■  7  f  ERAC  XA  +  4)  XB  +  (-h*  AC  +  E  BC  "  E  AB'  XA  XB 
+  3W(2)ac  (yA2  ♦  SA2)  ♦  31\’(2)Hr  (yR2  +  zB2)  +  30V(2)ac  +  wCz) 


-w(2V  (y.VB  + VbJ}  u.Pi(L)  '"Pi1 


...  (10) 


where  e\.  =  (8  -  fiwU)..} 


*  --  -y.  » *  w  '  *•  ‘  t  ^  g  \  ^ 


k  T 


xMGCA2/B2)  =  0.5-  [,|  CjCjBy  *  /(J  qCjEyj  2  -  4CACBCC  (E^) 


'  T  f  AC  *  EBC  -  EAB  ]2’} 

where  ij  =  AB,  AC,  BC  and  E-.  =  fswO)  -  6w(2)  ) 

J  ^  ij  ij  -B! 

k  TyBI,G(B2/L21)  .  3  ((!♦  Jv?  -  4W(2)ac  wfa)  K  -  [k«a(,  , 


...  (11) 


ij  WG 


'  JAb]  )b1VG  <CA  "  XA-*  tGB  "  'V  <CC  +  XA  +  XB'}  •  * 
whoro  «  .  OK*),^  (CB  .  CC  .  *A)  (CA  -  xA) 

*  (,S'(2)BC)  mCt  'CA  *  GC  *  XB>  <q»  -  xb) 

-  [w(a,AC  .  w(*>K  -  wtoj  (CB  -  x„)  (CA  -  xA); 


xind  xA  and  Xg  are  the  values  ol  the  order  parameter  in  nearest  neighbour 
positions  at  T  =  T  . 


5.  Calculation  of  T,  and  for  the  Fo-Al-Si  System 

The  following  procedure  was  adopted  to  obtain  optimum  values  of  and  W*-2) 

(BCC  Al-Si ) .  examination  of  various  pairs  of  solutes  shows  that  for  BCC  alloys  w^/W^s0.5, 
as  is  indeed  the  ease  for  the  other  pairs  concerned  in  the  present  case  (Table  I).  It  was 

therefore  assumed  that  N  "‘/K  ^  BCC  (Al-Si)  =  0.5.  Inspection  of  the  superlattice  intensities 
fifi*  (5)  show  that  a  considerable  degree  of  order  exists  in  the  B?  lattice  prior  to  DO 
formation;  the  maximum  theoretical  values  of  x.  and  xH  can  therefore  be  used  in  preliminary 
calculations.  .  A  B 

Trial  values  of  K  ,\i-si ( l^CC )  were  then  inserted  into  equations  until  a  reasonable 
fit  was  obtained  for  both  T^  and  T  for  a  selected  alloy  composition. 


As  shown  in  (Fig.  6)  there  is  a  limit  beyond  which  further  changes  in  N^(and  K^} 


are  counterproductive,  which  reflects  the  assumption  that  the  real  values  of  and  Xg  are 


lower  than  the  maximum  values  assumed  in  this  preliminary  calculation.  Varying  x^  and  Xg 

in  an  iterative  fashion  can  then  be  used  to  establish  more  realistic  values  of  x.  and  xD 

A  B 

as  shown  in  (Table  III).  Hus  result  can  be  cross-checked  for  other  trial  compositions  and 
a  pattern  established  for  the  variation  of  and  Xg  with  composition.  This  appears  to  be  a 

very  simple  variation,  with  x^  constant  at  0.18  (compared  to  the  initially  assumed 

maximum  value  of  0.25),  and  equal  to  a  constant  fraction  of  C.,  (the  constant  being  the 


l  equal  to  a  constant  fraction  of 

ratio  (xA/ x^(max) )  or  0.72).  These  values  of  ,  x.  and  x„  when  inserted  into 


A  B 

equations  (11  fj  12)  yield  the  results  shown  in  (Table  IV),  and  a  more  direct  comparison  with 
experimental  results  is  given  in  (Table  V),  and  (Fig.  7). 
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Trial  Values  ofW^-f  A  1-Si: 


*fefg.  6  Effect  of  varying  W'*1^  Al-Si  on  calculated 
i'x  and  ID  •t^iSpera'tui,QS-u§-i‘ng  tlie  a.ssuinption  of  full  A'2/|2  .ordering.  bt  T^>. 

f>.  Discussion 


The  degree  of  agreement  obtained  is  very  satisfactory,-  and  further  calculations  dir 
this- 'System  based  on  these  interaction'  paraiiietefs  should  yield  equal  ly  reliable'  results.- 
There* 3s;  some  indicat  ion  that  the  equi l  ibri uiii  values  of  and  x^.  for  <t lie  bd nan?  al  loys  Rc-Ai- 

Wi&fe^Si!  are'  in  fact  slightly  higher*  than  would; Do*  indicated  by  the.  simplifying  assumptjbitSK 

just  outlined  (see  Tabie  ii)  .  Whether  it  is  necessary  to  calculate  equilibrium  values  of  x . . 

”4*  ...  ,  i  i  ■*  *  '  **  J  1  ^ 

and' -jtg^gordusly  for  every  composition  instead;  of  Using  ;the  simpilfying; formula' depends  oir 

f  the  accuracy  with  which  the  ordering  temperatures  need  to  be  known  and,  also  to  some  extent  tons 
j^'ie  c6bl‘ing  rates  required  to  attain  equilibrium  under  .practical  conditions.  Table  A'l  gives, 
some  ‘idea* of  the  effect  of  small  variations  in  Ordering  parameters.  Additionally,  the  values;. 

of  T  and  T  shown  in  (Rig.  7)  refer  to  the  so-called  constitutive  phase  boundaries  and,; 

*X;  j,  * 


^analogous  to  the  Tq  lines  of  more  common  phase  diagram  calculations,  assume  no  two  pha 

regions  exist.  To  establish  the  full  phase  diagram-  a  more  complex  minimization  program  must  b< 
C*  utilized  (1,  2,  8)  which  is  rather  lengthy  but  requires  no  further  new  input  data.  It 
shduidValso 'be -noted  that  all  calculations  have  been  made  on  the  basis  that  the  temperatures 
of  interest  are  <above  the  magnetic  transition  temperatures  in  the  system;  if  it  is  intended 
to  investigate  ordering  below  the  Curie  Temperature,  additional  terms  have  to  be  introduced 


into  equations  (10-^12) :i  In^either  circumstance  the  quoted  values  of  W.f*f  and  WW  -can  be 
used  to  derive  heats  of  formation  for  solid  Al-Si  BCG  alloys  which  can  be  used  to  refine  the 
i|^|ciiiat;ion:  of  the  liquidus  and  sol  idus  fine  in  these  systems;  fl;7.',  l8),v 
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Fig.  7  Comparison  of  calculated  and  experimental 
ordering  temperatures  in  the  section  Fe^Al-Fe^Si 
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TAB LB  I 


OPTIMIZED  VALUES  OF  AND  W2  (k  units) 


! 


SYSTEM 

W 

Hi 

KATSNEL'SON 

USED  IN 

RATIO 

m 

(13) 

THIS  WORK 

Vi 

FeAl 

W, 

1005 

1090 

1080 

I 

1063 

FeAl 

W2 

539 

583 

547 

(W„/W.) 

4  1 

0.54 

0.52 

0.50 

Al-Si 

W1 

-480 

mm 

Al-Si 

W2 

-250 

HI 

oywj) 

H 

Fo-Si 

mm 

1P  -  - 

2010 

2010 

Fo-Si 

V2 

1000 

1067 

1000 

(W2/Kj) 

0.50 

1000 

_ 

0.50 

TABLE  II 

DERIVATION  OF  A2/B2  ORDERING  PARAMETERS  FOR  JFo .Al 
SHOWING  THE:  LEE  EOT  OE  CHANGING  W(1)(Fc-Al) 


}  =  1080 

W(:)=  540 

Tx  =  101 2  K 

T 

_JL 

(4  AG/6  xA) 

.17 

789 

-326 

.18 

794 

.  19 

798 

-95 

.20 

SOI 

+  80 

.21 

804 

+  321 

•  L  U 

807 

+665 

1005  K1 

*>  \ 

•  r  540 

calc  T  =  9001 

X 

*A 

T 

V 

(6AG/5*a) 

,15 

"77 

-93 

.16 

783 

-15 

.1? 

789 

+81 

.18 

794 

+  204 

.19 

79S 

+360 

.20 

801 

+  560 
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TABLE  III 


ai 


II 


TYPICAL  TRIAL  CALCULATION  OF 
ORDER  PARAMETERS  AND  T  FOR  A  TERNARY  ALLOY 


CONTAINING  75Fe,  13.5AI  and  ll.SSi 


S\ 

h 

T 

_Z 

AG/a* 

6  AG/5 x 

169 

-.080 

1243 

+42 

-41 

171 

-.081 

1244 

+25 

-23 

[71 

-.081 

1244 

+  39 

+  13 

172 

-.081 

1245 

-  4 

-41 

L72 

-.082 

1245 

+  8 

-  4 

173 

-.082 

1245 

+21 

+  32 

[73 

-.082 

1246 

-23 

-23 

Input  Energy 

Parameters 

BWG  Values 

in  K  Units 

1005*, 

540,  -480, 

-250,  2010, 

1000 

C 

a  1 eulated  T 

=  1516 

;e  N 

^  FeAl  here 

is  not  the 

X 

value  Finally  adopted 

TABLE'  IV 


CALCULATED  VALUES  OF  T  AND  T  (°K) 

x  v 


m 

T 

X 

Si " 

T 

>' 

(■V 

.25 

1890 

.24 

1458 

- .  007 

.  23 

1855 

.22 

1433 

-.022 

.21 

1816 

.20 

1405 

- .  036 

.19 

1773 

.18 

1370 

-.050 

.  17 

1724 

.  16 

1334 

« 

o 

O' 

in 

.15 

1671 

.  14 

1293 

-.079 

.13 

1613 

.12 

1247 

-.094 

.11 

1548 

.  10 

1196 

-.108 

.09 

1477 

.08 

1139 

-.123 

.07 

1397 

.06 

1074 

-.137 

.05 

1307 

.04 

998 

-.151 

.03 

1204 

.02 

908 

-.166 

.01 

1083 

0 

794 

-.180 

1080 

K 

540 

W 

18  (at  Ty) 


C„(-0. 72) 


t/MJfji.  V 


COMPARISON  OF  CALCULATED  AND  EXPERIMENTAL  VALUES 


T  (OK) 


T  (OK) 


•5i(at%) 


Exptl 


Calcd 


Exptl 


Calcd 


0.23% 

0.21 

0.195 

0.17 

0.15 

0.135 

0.115 

0.105 

0.085 

0.065 

0.045 

0.02 

0 


x.  =  . 18  x„  =  - . 72x  CD 
ABB 


(Calculated  value  for  T  Fe3Si  is  1890°K) 


TABLE  VI 


(a)  The  Effect  of  Arbitrarily  Reducing  the  Order  Parameters 
on  the  Calculated  T  Temperature  for  Fe75,  A1  13.5,  Si  11.5 


calc  T  (°K) 
v  v  / 


(b)  The  Effect  of  varying  the  order  parameters 
within  the  limits  found  for  other  compositions 


calc  T  ( °K) 

.  .  y 

175 

-.082 

1245 

18 

-.097 

1235 

16 

-.11 

1181 

19 

-.13 

1202 

18 

-.12 

1201 

aMSMiittMMI 


CALCULATION  OF  THE  FCC  +  BCC  PHASE 
FIELDS  IN  DISORDERED  Fe-A£-Ni,  Fe-A£-Mn  AND 
Fe-Ni-Mn  ALLOYS  BETWEEN  9006C  and  1200°C 

The  ManLabs  Data  Bank  was  employed  to  compute  the  fcc/bcc  two 
phase  fields  in  the  iron-aluminum-nickel,  iron-aluminum-manganese  and 
iron-nickel-manganese  systems  between  900°C  and  1200°C  for  disordered 
alloys.  The  results  are  shown  in  Figures  1-3.  Background  information 
and  details  of  the  calculation  are  presented  in  Tables  1-8.  The  calcu¬ 
lation  of  ternary  equilibrium  is  done  by  employing  CALPHAD  descriptions 
of  the  binary  systems  to  synthesize  the  ternary  fee  and  bcc  disordered 
phases  as  illustrated  in  Tables  1-8.  The  procedure  and  numerical  values 
have  been  published  as  referenced  below. 

1.  L.  Kaufman,  CALPHAD  l  7  (1977) 

2.  L.  Kaufman,  CALPHAD  2  117  (1978) 

3.  D.  Dew-Hughes  and  L.  Kaufman,  CALPHAD  3  175  (1979) 

4.  L.  Kaufman  and  H,  Nesor,  CALPHAD  2  325  (1978) 


TABLE  1 


LATTICE  STABILITY  VALUES  FOR  THE  ELEMENTS 
(Units  of  J/*ol  and  J/aolK  will  be  used  throughout) 


Elenent 


(L-Liquid,  P-Priaitive  Cubic-BMn.  K-Coaplex  Cubie-oMn) 


*Gb  -  *GPc‘ 


14744.416-8. 196ST 


*GL  -  *GbCC  «  13807. 2-7. 6316T 

*G,>CC  -  *Gfcc  ■  -5235*9. 4T-0.529S(E-2)T2*0.922t(E-6)TS 
*Gfcc  -  *Gbcc  ■  6108. 64-3. 4618T-0.7472(E-2)T2*0.512S4(E-5)TS 
*G^CC  -  *Gbcc  ■  S4S1. 752*7, 447S(E-S)T2-12(E-5)TS*20.5434(B-B)T4 

•GfCC  -  ‘Ghcp  .  1828.4  -  4.686T 


4.686T 


*Gfcc  -  *GP  »  -1151  -  0.837T 
*GP  -  'G*  .  -418  ♦  0.167T 

*GL  -  *Gbcc  -  16233.92  -  8.368T 
*Ghcp  -  *GbcC  .  -4351  ♦  3.7656T 
*Ghcp  -  *Gfcc  .  -3347.2 
*Gfcc  -  *GP  .  -4184 
•CP  -  *GK  .  2929 


Teaperature 
Rang*  K 

16«5<T 

1665<T 

1100<T<1800 

300k.T<1100 

(kT<300 

30CKT 

300<T 

300<T 

300<T<2400 

300£T<2400 

30(KT<2400 

300<T 


*Gb  -  *Gbcc  m  14644  -  9.623T 
•G1,  -  *Gbcp  .  9205  -  7.U3T 
*CL  -  *Gfcc  .  16401  -  10.878T 
’Gbcc  .  •CP  .  3975  -  2.887T 
*GP  -  *GK  .  2259  -  2.259T 

’  *GpCC  *  1 4 77 -0.5 HOOT- 2. 7420 (E-2)T2* l  ,6534(E-6)T3 
*G  CC  -  *G  •  611  ♦  13.101T-2.1240(E-2)T2*0.8J96(E-5)T3 


•GL  -  *GfcC  .  10711.04  -  U.S06T 
*Ghcp  .  *Gbcc  <  -4602.4  ♦  3.01248T 
*Gbtp  .  *Gfcc  .  5481.04  -  1.79912T 


*Gbcc  -  *GP  .  -837 
»GP  -  •G*  »  -4187 


1220<T 

1220<T 

1220fT 

122IKT 

12205.T 

400<T<1220 

400<T<1220 


•GL  -  *Gfcc  «  17614.64  -  10.209  T 
*GbCC  -  *GfCC  •  5564.72  -  1.046  T 

*c*CC  -  *GfCC  •  3932.96*4.1087  E-3  T*-0.48S3  E-S  TUMI  E-9  T1* 
•GhCp  -  *GfcC  •  1046  ♦  1.2552  T 


1000  <T 
0  -  1800 


10-  «  Ea 


76 


TABLE  2 

THE  IRON-MANGANESE  SYSTEM 


ANALYTICAL  DESCRIPTION  OF  THE  IRON-MANGANESE  SYSTEM 


¥ .  Vi, 

*FaW'1M?4*Ft'2,M9W 

\.W*4,0°,F.‘S1MW 

.^(-UITO) 

-1.4602x10’ *7*) 

«f.*»I-wm»f.:sm‘Sr,I 


V  .  S4  «F,  S4.  Vs*  Co^xiUon 

,*t,F»t,uW,uW 

-xf,V»6.M?  »<**<• 

X^I-A.EIEX^-J.SOE^J  04*^4  1 


-,f,^IA.M7 


.^[-n.Htx^S.EITx^  04X^4  1 
•>.*>04xl0’*Tl 

W-7-,,2*f.’2*l#W  #*V*‘ 


1S004T4  JOOO 
’Reftri  to  liquid 

1004  T  4  >000 
'Rofori  to  bee 

1004  T<  >000 
'Roftrt  to  fee 

*004  T4  1100 
'Roftra  to  Ma 


4004  T4  1>00 

'Rofori  to  4h 


EXPERIMENTAL  THERMOCHEMICAL  DATA  FOR  I  RON -MANGANESE 


*Mn 

EHfcc 

ALLOYS  (t) 
Ejjfec 

ECL 

0.1 

* • • • • (14S0X) -  -  ■ 
•1314 

•i.2i* 

ll«6lX) 

397 

0.2 

-2456 

•2.297 

703 

0.1 

-1431 

•3.167 

924 

0.4 

•4204 

-3. *37 

1054 

0.$ 

•472*11674 

•4.2S$iO.I37 

110011674 

0.6 

•492$ 

•4.377 

1054 

0.7 

-4690 

•4.113 

924 

0.1 

•1*95 

•1  111 

703 

0.0 

•2319 

•2.054 

197 

(a)  Calculated 


TABLE  3 

THE  IRON-ALUMINUM  SYSTEM 


V.s*-«p/s*#.«ArsJ|  0im,nu 


IX...42760.1..94212)  l«.J. A02-X..11.014]  0<x..«1  400<T<1900 

Ft  Al  H  Al  M  Al  ft  At  At  ‘MfoM  t»  UquM 


•WA|l't,».Ml™|  ****!**  ^iwbcc 


WMim'\lUnt*'  WA>tV4I“-*M,'-,I*>  "“At*1 


Coapoiltlon  Cium 


P.2I0  M0.710 


•1114*  -I, 111  iJj-0.710 


stiMuiv  or  (xmiWKTU  twunocwmical  mta  k*  imn-auminun  auo»s 

(1-x)94(U«iM)  •  xAl(ltqul4)-»  Nf,  .Al  (tHul<l)(tl7ir) 

AC  AM  >  AC  I  AO 


0.) 

•90401 2100 

•411012000 

0.4 

-2144012100 

0.7 

-2001012100 

0.2 

•14910 

•moo 

0.1 

-22111 

0.9 

-14110 

0.3 

•19910 

•14020 

0  4 

-22120 

0.9 

•MM 

(!•«)  F»WC  *  «  A1,m  ♦ 


*  AC 

4H(29U) 

t 

4G(U71KJ 

4H|290() 

1  AC[II7») 

0.1  -977011100 

•177912000 

0.4 

•2470012100 

-2000012000 

0.704  -2490012100 

0.2  -14400 

•11000 

O.S 

•21910 

-21100 

0.710  -24110 

0.1  -21700 

•11710 

i.  Ml 

•21410 

•24110 

■tUNllM 


U4Ul*l  ^L) 


f4>- 


(a)  larutatcd 


Ibf  Ob>erved 


/ 


TABLE  4 


ANALYTICAL  DESCRIPTION  OF  THE  IRON-NICKEL  SYSTEM 


Phase 

*"*  ■  %  '  *F.  •"?.  •  *Ni  ‘"Hi 

S  S  xFe  SFe  xNi  5Ni 

-R[xpa«nxFe*xHilnxN4] 

Composition 

Range 

Comments 

Liquid 

xPexNi  t*8368  Xpe*32217  x^j) 

xFexNi^’2,72xFe*9,293xNp 

0<  xNi  <X 

0  <T  <18001 

fee 

xFexNilxFe(2092+3,*314xl0’*T* 

-3.2652xlO**T’)+xNi(-34881 

*2.4404x10*2T  ♦  2ii802xl0’s7,)J 

xFexNi^xFe(7,6628xl0’3T'4,8990 
xlO‘*T*)  ♦xNi(-4.8808xl0',T 

♦3.1203x10- *T2) } 

0<  xKi  <1 

0  <T  <1800K 

bee 

xF3XNi^xFe^339*3,327Sxl0",Tl 

♦3.1740xl0-*T‘)*xNi(-l6276 

-3. 5412xl0*,T**2. 6382x10' »T*) ] 

xFexNi  lxFe^'2,655x10' *T*4*76l 
xl0‘*T! )  ♦xNi(-7.0824xl0-*T 

♦  3. 9573x10*  *T* ) ) 

0<  XNi<  1 

0<  T  <1800 

Compound 

*"  •  «  »?.  •«£'  •  »s,  *hJ;c 

iS  .  S  •  x*„ 

Composition 

Comments 

F«b.25Ni0.75 

AuCUj 

Structure 


-15592 


-10.83 


xNi 


0.75 


FIG.  38  (observed)  (40) 


FIG.  (calculated)  (30) 
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TABLE  5 

ANALYTICAL  DESCRIPTION  OF  THE  ALUMINUM-NICKEL  SYSTEM 


>H*N 


Composition 

Range 


Liquid 

j  xAlXNi 1 XHi  <"203376“  1  •  3924*  •  0~*  TJ  +5.4626*  1 0'5T*; 
*XA1 ('9208 1 ' ' ' 3924* 1 ^  T*  4  5 ’ 4 626x 1 °" * T ’ 3  3 

XAlXNi  {xNi ( 1 38 • 1 5-27  * 848  *  °"iT*d  *  • 939  *  °"*fl 
+xA1  (130. 206-27 .848  *1 0~*T+8 1 . 939  *1  o"*T* )  ] 

0  <*„.  4. 

fee 

xAlxHi  lxNi  1 89987“ 1  •  392h“!Q_1Tj  +5 .4626*  lO*5?* 
+xA1 (- 1 9907- 1 . 3924“ IO-1T,+5. 4626* 1 0"*T* ) ] 

xAlxNi^xNi^  *49.452-27.848><l0"’lT+8l  .939xlo"<I8) 

♦XA1 ( 182 . 506-27 . 848“ 10_J T+8 1 . 939* 10"*T* ) 1 

0  <x  <, 

Ni  1 

bee 

XA1XNIS  xNi  ("220984" 1  * 3924*  •0~,T**5 . 4626*  JO-*!1) 
♦xA j (-28 1 662+0. 92048* 1 0"‘ T* - 3 . 6 1 1 6* |0_,T 1 ) ] 

*AlXNi f  XNi  < 14 1  •  °84-27 . 848*  1 0_lT+8 1 . 939“  1 OV ) 
♦xAl (-220. 622+18.410*10'*  T-54 . 1 74* I0"*T* ) ] 

0  <xNi  <1 

Compound 

AWxAl‘HAl'xNi%i 

AS'S'XAl‘SArXNi*SNi 

Composition 

Al0.75Nl0.25 

-33465-2.611  x  10'V  ♦  1.0242  x  10'V 

20.00  -  5.222  «  10*2T  ♦  1.5364  x  10*V 

x"i*°“2S 

A10.6Ni0.4 

-49352  -  0.33418  x  10'V  ♦  1.3110  x  10*STS 

26.84  -  6.683  x  10*2T  ♦  1.9665  x  10'ST2 

x^-0.40 

Al0.25Nl0.75 

|  -38878  -  2.611  x  10* V  ♦  1.0242  x  10*V 

1  , 

27.57  -  5.222  x  10*2T  +  1.5364  x  10’V 

J*  n 

C 

© 

•“J 

Consents 

300  <iT  42000K 
defers  to  liquid 

300  <T  f 2000  K 
•Refers  to  fee 

300  S<T  £2000  K 
'Refers  to  bee 

Comments 

•Refer*  to  fee 

•Refer*  to  fee 

•Refers  to  fee 


If  ! 


TABLE  6 

ANALYTICAL  DESCRIPTION  OF  THE  MANGANESE-NICKEL  SYSTEM 


Phase 

* 

_  it$  v  Ini  v  Om$ 

H  •  “m  xMn  “Mn'^i  “ni 

ES*  =  S4-Xj.  °S*  -x„.eS4 
^4n  Mn  Ni  Ni 

+RtWnxMn+xNi*nxNi] 

Composition 

Range 

Comments 

Liquid 

-xMnXNi&4434xMn+76986xNi] 

-WNi10’878 

0<xNi<l 

1200<  T ■«  1800 

0  Refers  to  Liquid 

bcc 

-xMnxNit4S606xMn+S31S8xNi] 

■XMnXNi3-640 

0<xNi<l 

1200  <  T<  1500 

0  Refers  to  bcc 

fcc 

-xMnxNit51882xMn+64434xNi^ 

'WNi10-878 

0<xNi*1 

300  <  T<  1800 

0  Refers  to  fcc 

P(BMn) 

-xMnXN#3178xMn+39330xNil 

'xMnxNi17,991 

0  *  xNi  *  1 

700  <  T  <  1500 
°  Refers  to  gMn 

K(oMn) 

•xMnxN#3178W39330xNil 

'xMnXNi17,991 

0.<xNi.<l 

400  <  T  <  1200 

0  Refers  to  aMn 

Compound 

AH  ‘  H-xMn°llSn-xNi0,,!i 

iS  *  8"xMn°SMn’xNi°SNi 

Composition 

Comments 

Mn0.750Ni0.2S0 

-27379 

-17.730 

x*.  =  0.250 

0  =  K(aMn) 

Mn0.667Ni0.333 

-14028 

+1.302 

xNi  "  °’333 

0  *•  fcc 

Ktn0.500Ni0.500 

-18891 

♦1.360 

x+A  =  0.500 

0  •  bcc 

m 


Mn0,333Ni0.667  '16777  -°-557  x5i  =  °'667  9  '  fcc 


MnQ  2soNio  750  -15708  -’.040  x*.  =  0.750  0  =  fcc 


TABLE  7 

ANALYTICAL  DESCRIPTION  OF  THE  ALUMINUM-MANGANESE  SYSTEM 


EH*<VHArxMn°<  Es*=SSrSArxMn‘SMn 

+R[xAl£nXAl+XMn£nxMn] 


Liquid 


-xAlxMn57521 


'xAlXMn 


Composition 

Range 


0<x..  <1 
Mn 


0<x..  <1 
Mn 


800<T<1800 
“Refers  to  Liquid 


1000<T<1700 
“Refers  to  bcc 


TABLE  8 

DESCRIPTION  OF  THE  FREE  ENERGY  OF  TERNARY  PHASES 

Definition  of  the  free  energy  of  the  bcc  phase  of  the  Fe-A£-Ni 
system  x*x^,  ysyNi  can  readily  be  stated  as 

Gbcc  =  ( 1  -x-y)Gb^c+xG^c+yGb^c+RT  [  ( 1  -x-y)£n  ( 1  -x-y)  +x£nx+y.£ny  ] 

+  [ ( 1  -x-y) BFEAL+xBALFE] ♦  y-  [(l-x-y)BFENI+yBNIFE) 

+  -^-[xBALNI+yBNIALj  (1) 

x+y 

Similarly  the  free  enorgy  of  the  fee  phase  can  be  written  as 
Gfcc  *  (l-x-y)Gp^c+xG^c+yPjjj?c+RT[(l-x-y)£n(l-x-y)+x-6nx-*y£nyj 

+  [  ( 1  -x-y) AFEAL+xAALFH]  +XililXlZ.  [(1-x-y) AFENI+yANIFE] 

>  [xAALNI+yANIAL]  (2) 

x+y 

On  the  basis  of  Eqs  (1)  and  (2)  the  partial  molar  free  energies 
of  iron,  nickel  and  aluminum  in  the  fee  phase  and  the  bcc  phase  can  be 
written  for  the  orientation  Fe-A£-Ni  considered  above.  Thus  with  x^x^ 
and  y=yNj  the  partial  free  energies  in  the  fee  phase  are  given  by 
Eqs  (3) -(5) 

GPeC  =  GR}C  +  R^nO-x-y)  +  x2AFEAL  +  y2AFENI 

+xy  (AFEAL+AFEN I  -  AALN I )  +  x2[|j~^-  -  ]1  (AALFE-AFEAL) 

+y2[0^  '  11  (ANIFE-AFENI)  -  (ANIAL-AALNI)  (3) 

GA£C  =  GaT  +  RT^nX  +  0“X)2AFEAL  +  y2AFENI 

-y ( 1  -x)  (AFEAL+AFENI  -AALN  I )  ♦  " 1 ] x2  (AALFE-AFEAL) 

2 

♦y2(y^-  -1]  (ANIFE-AFENI)  +  [y(l-x)-x2XANIAL-AALNI)  (4) 
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and 

SNiC  =  GNiC  +  Rltny  +  x2afeal  +  Cl-y)2AFENI 

-x(l-y) (AFEAL+AFENI-AALNI)  +  x2[~ --  1] (AALFE-AFEAL) 

+y2[O^F  +  ^  (ANirE-APEND 

+  Ix^  t2x+y'Xy'y2l(ANIAL_AALNI)  (5) 

Similar  equations  can  be  written  for  the  partial  free  energies  for  iron, 
aluminum  and  nickel  in  the  bcc  phase  merely  by  replacing  the  symbols  fee 
and  A  by  bcc  and  B  in  Equations  (3) -(5). 

Tables  2-7  contain  the  explicit  equations  for  the  A  and  B  values 
in  all  of  the  six  binary  systems  of  interest  to  compute  the  fcc/bcc 
equilibria  in  Fe-A£-Mn,  Fe-Ni-Mn  and  Fe-Mn-Af.  In  addition,  Table  1 
contains  the  free  energy  differences  between  the  fee  and  bcc  forms  of 
Fe,  M,  Ni  and  Mn  so  that  the  relative  free  energies  of  the  fee  and  bcc 
phases  in  each  ternary  can  be  specified.  Simultaneous  equilibration  of 
the  partial  free  energies  of  each  component  in  the  fee  and  bcc  phase  permits 
calculation  of  the  tie  lines  across  the  two  phase  fields.  The  latter  are 
shown  in  Figures  1-3. 
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